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SUMMARY 


Sound wave propagation in a soft-walled rectangular duct with steady uni- 
form flow was investigated at exhaust conditions. An analysis was developed 
incorporating the solution equations for sound wave propagation in a rectan- 
gular duct x^ith multiple longitudinal x^all treatment segments and uniform 
flow* Modal analysis x-iras employed to find the solution equations. The 
analysis was used to study the effectiveness of a uniform (one-segment) and of 
two-sectional liner in attenuating sound pox^er in a treated rectangular duct 
xi^ithout flow (M = 0) and x^ith uniform flow of Mach 0.3. The optimization 
studies were, made for a duct height- to-x^ave length ratio of 1*535 and x^all 
treatment length- to-duct height ratio of 3.43. The optimum predicted sound 
attenuation xj as . compared x^ith measured laboratory results. 

Several two-element treatment configurations were designed to reproduce 
the optimization conditions and then tested. Good correlation x*7as obtained 
between the measured and predicted suppressions x^hen practical variations in 
the modal content and impedance were taken into account. Maximum suppression 
was found to be very sensitive to the impedance components of the liners: a 
variation of any one of the normalized resistance or reactance components by 
4- 0.1 can lox^er the predicted peak suppression of 35 dli by 5 to 10 dB. 

The analysis was then employed in a series of parametric optimization 
studies, first, of the effect of the ratio of duct heigh t-to-x^avelength, n, 
on the optimum sound suppression in an acoustically lined rectangular duct 
and, second, of the effect of variation of liner segment length ratio on optimum 
sound suppression when the total length of the two-element treatment is kept 
constant. T?^o-segment optimized liners wave shoxi/n theoretically to increase 
the attenuation of sound as compared to an optimized uniform liner. That 
suppression increase depends upon the r\ value, i.e., upon the frequency, and 
incident source mode distribution. The most practical use of this enhancement 
occurs over an approximate range of frequency parameter n (H/X) from 1 to 5. 
Maximum sound attenuation due to a tx-jo-segment liner of a given length depends 
upon the relative length of the liner segments. Optimum suppression increases 
x^ith increase from zero to slightly over one in the ratio of thp length of the 
upstream liner segment to that of the dox*7nstreara segment. Maximum optimized 
suppression is attained at, or slightly above, the segment length ratio of 
1.0; i.e., when the two liner segments are approximately of the same length. 



SECTION 1.0 


INTRODUCTION 


Substantial noise reductions for modern turbofan engines have been 
accomplished with the use of acoustic treatment, or sound absorbent liners, on 
the walls of the engine Inlet and exhaust duets.; The state of the art of 
acoustic treatment design has developed to the point where, if further gains 
are to be achieved, it is necessary to investigate detailed flow duct acoustic 
phenomena which cause one engine configuration to differ from another. These 
effects include the nature of noise generation mechanisms, acoustic propa- 
gation through flowing air, and the optimization of liner sections, among 
others. Such phenomena are not susceptible to the empirical approaches which 
have served to develop past designs* An excellent survey of current jet engine 
acoustics is presented by Nay f eh, Kaiser, and Telionis^^X* 

The basic theory of acoustic propagation in ducts \^ith airflow has been 
known for sometime, but the correlation of the theory with experiment has been 
scarce. The emphasis of this study is on the application of modal analysis 
of wave propagation in ducts, with uniform airflow, to design acoustic liners 
which are optimized for maximum suppression under a given set of constraints* 

A design method is developed and the theory is tested by comparing analyti- 
cally predicted suppressions with measurements made in a laboratory duct* 

The model for the theory and experiment was chosen to he a rectangular 
duct* This simplifies the analysis and allows experimentation under control- 
led conditions* It is not without practical application, as well, since 
propagation in a hi^ radius ratio annulus is approximated closely by propa*^ 
gation in a rectangular duct(^^. A major objective of the investigation is to 
deteinhine how detailed the analysis must be to predict accurately . the physical 
phenomena. 


1.1 PROGRAM OBJECTIVES 

The primary objective of the program was the analytical and experimental 
evaluation of the multielement suppressor concept* The evaluations were 
conducted at conditions of aeroacous tic flow fields representative of engine 
fan exhaust ducts. The program was accomplished in two major Tasks, Task I 
being an analytical study and Task II an experimental program* 


Superscript numbers refer to references contained in Section 7*0* 
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1.2 TASK X OBJECTIVES MD SCOPE 

The objective of the analytical study in Task I was to develop a method 
for finding the optimum design of a two-pane l-element exhaust durt suppressor 
employing step discontinuities in the wall impedance, given a set of specific^ 
but arbitrary, Initial conditions* The method makes use of the theory of 
modal propagation A^ath airflow in an acoustically treated rectangular duct; 
specific conditions for the analytical method are: (1) the duct is two- 

dimensional, (2) sound propagates with airflow of uniform velocity across the 
duct, and f'") the sound source consists of a discrete frequency, statistically 
stationary ‘Tith time* A further objective was to use the resulting analysis 
to determine an optimum liner configuration under a specified set -f conditions 
to build and test the design; in a laboratory duct (under Task Tl) , and tb 
compare predicted with the measured suppression results. Finally, the scope of 
this task included using the analysis in two parametric studies to deterinlne 
the amount of suppression that can be obtained with two panels having optimized 
impedance components, one to investigate the e£fe:ts of the due t- height /wave- 
length ratio, and the other to determine the effect of the ratio of the 
lengths of the two liner sections. 


1*3 TASK II OBJECTIVES AND SCOPE 

The objective of the experimental program in Task 11 was to demonstrate 
experimentally the optitnized treatment design determined in Task I and provide 
data for correlation with the analysis* The test program was accomplished in 
two main phases, with a short preliminary test preceding both ol these. 

Prior to the initiation of the parametric design optimization in Task I, 
it as necessary to measure the modal source characteristics of the duct 
under the optimization conditions. This was accomplished in a preliminary 
dtict test in a hardw all configuration* 

Upon completion of the design optimization in Task I, the Phase 1 test 
program was initiated. Six treatment panel hardware designs were determined 
which \would .■■bracket” the desired impedance components. These configurations 
were then tested in the duct for third octave transmission loss over a wide 
range of frequencies and for narrowb and transmission loss at the design fre- 
quency. The treatment configuration showing the best suppression was then 
subjected to detailed modal, measurements. 

The objective of the Phase II test program was to improve upon the results 
of the Phase I testing, based on evaluation of the Phase I test results. This 
led to ah Analytical reopt iMzation^ b on revised source modal character-- 

istics and test of the revised design. In addition, such effects as sensiti- 
vity to source modal characteristics and wall impedance were investigated in 
detail as part of Phase II. 


ORIGINAL PAGE IS 
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SECTION 2*0 


ANALYSIS 


Starting with the wave equation in Cartesian coordinates, modal analysis 
is employed to find the solution for the wave propagation in a sof t'-walled, 
multisectional rectangular diict with flow* 



The sound field in a homogeneous medium with uniform steady flow is 
described by the following linearized wave equation (3). 


,2 1 j. , m2 



where ■ 

the Laplacian operator 
pCxjZjt) - the acoustic pressure 
c = the speed of sound 

K - Uniform steady flow Mach number in the positive axial 

direction, z 

Two-dimensional: duct, g^metry is assumed, i.e*, the rectangular duct is 
modelled in two dimensions as shoim in Figures 1 and 2 for a single acoustic . 
treatment segment* The sound pressure waves enter the treated portion of the 
duct at z = 0 and propagate in the positive z-direction* A positive flow ^ 
Mach number in steady uniform flow in the positive z-direGtion and 

thus simulates exhaust conditions, A negative Mach number xndiGates flow 
opposite to that of the propagatlGn of sound, i* e* , simulates inlet condxtions . 

V The eliminatiph of time dependence , /VTbich is assumed tp be of the. form . 
e-itiit ^ yields the two-dimensional , time- independent wave equation , 


where 




X and z are Cartesian coordinates (See Figure 1) 






Plane j Plane j 1 



Figure' 2 . Geometry ol Duct Section. 


k = — is the wave number in free space 
c 

w “ 2frf is the angular frequency 
f = frequency, Hz, 

The solution of Equation 2 for the m-th mode of the acoustic pressure in 
a single uniform duct section is given by the sum of the incident and reflec- 
ted waves for that mode, i*e., 


. J 


o (y+ 


m j+1 


The first term in the brackets represents the forward traveling (incident) 
wave, while the second one is the backward traveling (reflected) wave. The 
+ sign on the superscripts refers to the former, the - sign is associated 
with the latter. 

The spatial part of the complete solution for the acoustic pressure is 
the sum of the characteristic solutions 


iK^(z-z . ) 




iK (z-z ) _ 

e " ^ f) 

Til m ti 


where 

X = the transverse distance 

Y|jj = transverse mode eigenvalue 

m = transverse mode number 

NM = number of transverse modes in the expansion 

j, j+1 = axial station indices 

H = duct height 

= modal expansion coefficient 
= eigenfunction for the m-th mode 

± = superscipt signs denoting the direction of wave propagation 

The axial propagation constants k for the forward and backward propa- 
gating mode m are related to the flow Mach number M and transverse eigenvalues 




k 


+ 

*■ 

K e 
HI 



( 5 ) 


The normalized symmetric modes are given by 


4i (x) = 


(l+cos y ) cos(y ^ ) + sin Y sin(V 

ta m ri tn m H 


( 6 ) 


H( l+cos Yjj,) ^1+ ™ 


siny 


in 


m 


and the normalized antisymmetric modes are given by 


i|> Cx) = 
m 


-(1 -cos y^) cos (v^ ^ ) + sin y^ sin(y^ |) 




H(1-cosy ) ( 1“ 

'm V — - 

' m 


(7) 


The expressions in the denominators of Equation 6 and 7 are the normal- 
ization factors obtained from the integral 


CNORMj„ = 



(x) dx 


( 8 ) 


For symmetric modes, it yields 


CNORMj^ = yH(l + cos}^,)(l + 


(9) 
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and for antisymmetric modes 


CNORM^ = yiKl - cosymHl - Hiny^i/yni) 


( 10 ) 


Th3 eigenvalues and the eigenvectors are determined from the wall 
boundary conditions which can be expressed in terms of impedance ratios or 
admittance ratios* 


The boundary condition to be satisfied by the symmetric modes is 


/3kH = 



( 11 ) 


and by the antisymmetric modes is 


PkH 



± 

1 + cosy 

m 


sin Y 
m 


( 12 ) 


where , 

6 is the duct wall admittance ratio (based on the e time periodicity), 

n = exponent depending on the type of boundary condition 
We set n = 1 for the particle velocity continuity wall boundary condition 
and n - 2 for particle displacement continuity wpll boundary condition. 


2,2 RECTANGULAR DUCT EIGENVALUE DETERMINATION 


The determination of the transverse eigenvalues requires finding the 
roots of two complex transcendental Equations ^ (11) and (12), This can be 
done by graphical or analytical methods. The mapping of the variable 6kH 
in the complex eigenvalue plane for symmetric and antisymmetric modes, at 
Mach number M = 0.0, is sho^m in Figures 3 and 4, respectively. 

For every (physically possible) value of 31di there is a corresponding 
sequence of symmetric and antisymmetric eigenvalues. Each eigenvalue 
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Real {>') 












in the sequence falls within one of the regions delineated by the cross- 
hatched lines separating the tuode orders. While for Mach 0.0 the map applies 
to all kH values, for Mach numbers different from zero a separate mapping is 
required for each Mach number and every value of kH. 

The roots of Equations (11) and (12) are determined analytically by 
employing an iteration procedure referred to as Bailey's iteration method^ . 

Given an estimate Yj of a root of the equation 

F(Yj) = 0 (13) 

Bailey's iterative formula for computing a refined approximation Yj+i of 
the root from the approximation yi is 










2F’(Yj) 


(14) 


where j is the iteration index, and the primes denote derivatives with 
respect to y* 

Bailey^s method is similar to the more familiar Newton-Elaphson iteration 
technique but it has a higher order of convergence, i*e*, it converges at a 
faster rate than the former. 

The success of an iteration scheme such as this is measured by its 
ability to rapidly converge to the correct eigenvalues* In the present case 
the eigenvalues are multivalued functions of the wall admittance and the 
convergence to the correct root is critically dependant upon the assumed 
initial value of the root at the start of the iteration. In the past, it has 
been customary to start the iteration from the sequence of hard wall roots on 
the real axis of the complex eigenvalue plane. Although these initial values 
start in the correct modal regions, there is no guarantee that the iteration 
path of convergence will not cross a branch cut separating modal regions. 

This may lead to double eigenvalues and/or missed eigenvalues, since the 
technique always converges to the eigenvalue nearest to the starting point, 
regardless of the branch cuts. 

It is Important that the eigenvalue routine be as accurate as possible 
in a program used as an optimization design tool, particularly when many 
cases are to be run in succession. A special effort was therefore made to 
develop an eigenvalue routine free of problems usually encountered in pro-’ 
grams of this type. 
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The adopted procedure consists of dividing the modal regions in the 
eigenvalue mapping plane into subregions, and choosing an initial value for 
each mode in the appropriate subregion* The subregion boundaries are either 
lines of constant magnitude of or lines of constant phase of 3kH. Typi- 

cal examples of the partitioning and initial values are shorn in Figures 5 
and 6 for the lower part of the symmetric and antisymmetric mode planes. 

Using these initial values, the iteration procedure is then performed at 
M = 0 until convergence criteria are met. For a nonzero Mach number, the 
Mach number is then increased from zero to its final value in small incre- 
ments, using the converged eigenvalues from the previous Mach number as 
initial values in the next step. This is done for both positive and negative 
Mach numbers at the same time, since the respective eigenvalues are needed 
for both forward and backward traveling waves. 

The assumption used in the above procedure is that the eigenvalues shifr 
continuously as a function of Mach number. Results of the program and previous 
experience in modal mapping appear to confirm this continuity even though the 
mapping of cases with the flow may be different from the M = 0 case, partic- 
ularly for low values of kH and continuity of particle displacement boundary 
condition. Checkout of the eigenvalue routine for several 3kH values and 
flow Mach numbers has indicated it to be a reliable, accurate, and rapidly 
convergent calculation. It is not infallible, however, and certain combina- 
tions of parameters may still cause convergence problems. Since the accuracy 
of the solution is very sensitive to the eigenvalues, each case must individ- 
ually be checked for convergence or the occurrence of double eigenvalues. 


2 . 3 DUCT SOURCE DETERMIUATIOM 


The duct propagation analysis requires knowledge of the pressure modal 
distribution over at least one plane in the duct. The complex modal coeffi- 
cients must be obtained, such that they contain both relative phase and 
magnitude information. 

If it is desired to compare the results of the program with a rectangular 
flow duct experiment, it is necessary to perform a modal measurement at the 
planes designated as source planes. This consists of measuring acoustic 
pressure cross-spectra between a stationary reference microphone and a micro- 
phone which can be traversed across the duct. The cross-spectral density, 
at a given frequency, as a function of distance across the duct, is then 
expanded in duct modes as 


S (to , x) 




|)+ aT <p. (y. 



(15) 


If the modal measurement is performed over a hardwall section of the duct, 
the forward and backward eigenvalues are the same, so that Equation (15) 
can be written 
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Figure 5, Initial Eigenvalues for 
Symmetric Modes. 


Real in 



Imaq (^) 


Figure 64 Ini"tia3. Eigenvulues for 
Antisymmetric Modes, 




(16) 


S(u,x) 



(At + Ap (jir f) 


Thus, the coefficients which present the sum of forward and backward tro.veling 
waves can be determined by standard Fourier methods. 

The system matrix can be written to allow this Interpretation of the 
source vector as a superposition of forward and backward traveling waves r 

The eigenvalues with flow over a softwall section of duct are different 
in the f onward and backward direction, so that in this case a variation of 
the procedure is required. We multiply Equation (15) successively by 


and integrate from zero to H, This gives two sets of equations. 


/ |)dx = ^ [a| (rt " ‘I /o%<n 1’^ 


and 


£)<!« . [a; (y[ P'fiCr; + A- X 


o •>•1 '>'1 


r X 


^) 4 > (/■ 


A A 


l>^] 


(19) 


vrhich can be solved for the A. and the A., 
out. J J 


since the integrals can be carried 


The theory behind the modal measurement requires the experimental deter- 
mination of the complex acoustic pressure profile as a function of x across 
the duct at given frequencies of interest. At a given frequency, the modal 
coefficients are determined by a modal decomposition of S(oi,x) in Equation 
(15) using standard Fourier analysis-type techniques. The magnitude and 
phase of S((iJ,x) are determined at a number of immersions by measuring the 
cross-spectral density between a moveable probe in the duct and a reference 
microphone fixed in the wall of the duct. The experimental setup for a 
rectangular duct is shown schematically in Figure 7, 
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Figure 7* In-Duct Modal Measui^ejnent Probe Configuration. 


The cross-spectral density of two signals is the product of the Fourier 
(time) transforms of the two signals in the form 


S 

rp 







(Xp,w) 


( 20 ) 


where 


Pp(Xp,oj) = 



reference microphone pressure signal 
probe pressure signal 

positions of reference microphone and probe 
angular frequency » radians /second 


The bar denotes Fourier transform with respect to time, and the asterisk 
denotes complex conjugate* If we make the modal measurement in a hardwall 
section of duct, the Fourier transform of Equation (16) gives 


p(x,tu) = 


C (O)) (y. f) 

j=0 J J J H 


( 21 ) 


where 


C. 

J 


-Aj + 


A. 

J 


Combining this with Equation (20) , we get 


( 22 ) 



foO I 

r 00 -1 

s (0;,x) =: 
rp. 

_n:=0 

'T' ^ 

>- 

4 

1 


(23) 


The first sum on the right hand side is constant with respect to x, so that, 
if we let 


B . (uj ) 
J 




(24) 
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we caa write 


CO 

S^^CUJ.x) ^ 2 Bj . ) (25) 


The left side of this expression is obtained from measureTnent, and the B, are 
then obtained from analysis . Since the Bj are proportional to the C ^ , tne 
relative magnitude and phase of the modal coefficients can be determined in 
this manner. • . . • 


2 A MATRIX SOLUTIQ]^ METHOD 

The two-dimensional characteristic waves propagating in the duct combine 
to yield the resultant x- 7 ave (Equation 4). The amplitude of the resultant 
wave is the complex sum of the amplitudes of the characteristic waves. It is 
convenient to represent multimodal acoustic propagation in a roultisectioned 
duct in the form of matrix equations. This is similar to the approach 
employed by Zorumski^^»^\ 

In the multielement duct, the finite duct sections are interconnected* 
Correspondingly, the characteristic duct solution (Equation 4) applied to 
each finite section must include acoustic coupling effects between the 
respective duct sections, i.e., the propagating wave must include trans- 
mission and reflection effects at . the appropriate, transverse planes of 
longitudinal wall impedance discontinuities. Figure 8 shows schematically 
a two-element duct with the indicated transverse planes. 

The condition for acoustic T^ressure continuity across a plane of axial 
impedance discontinuity at Zj = is, in vector form, 

{A'J} + {A+Jy 

(26) 

{^e+1) (4*(3+1),+ ^ 

For this a;xial compoTient of acoustic velocity, the continuity condition is 

.. ,, ( 27 ). 

2# Z 


18 



Plane Indices 



7 I Interface , 

‘ Plane 

“ ^j+l ^j+2 


Figure 8, Plane Designation Convention, 
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where, with the sign of the superscripts indicating direction of wave propa- 
gation, 

4 ' 4 

^ vector of eigenfunctions, plane j 

— "1 

{A ■^ } = modal solution vector at plane j 

T = superscript designating the transpose of a matrix. 

The elements of the axial impedance matrices are given by 





(28) 


Inhere 6^^ is the Kronecker delta. 

= 1, for 1 = J 
0, for I J 

By multiplying Equations (26) and (27) by 

} (29) 

and integrating across the duct, we obtain 

(A+«+l>> . [i-ifj+D.+J) + [r+«+1).-(3+D] {a-«+1>) ,3 


where [T] is a transtnission matrix from planes j to j+1 and [R] is a reflec- 
tion matrix at plane j+l. Thus, this equation relates the forward traveling 
wave at plane j+1 to the combination of the forward traveling wave transmitted 
through the discontinuity from plane j and the reflection of the discontinuity 
from the backward traveling wave at plane j+l. 

Similarly, by multiplying Equations (26) and (27) by 

(31) 

and integrating across the duct, we obtain 
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(32) 


,a- 3, , [T-i-tJ+D] (A-<J+1>) + [R-i+3] (A+J) 


Equation (32) relates the bacla^ard traveling wave at plane j to the backward 
traveling wave at plane j+1 transmitted through the discontinuity and the 
foni/ard traveling wave at plane j reflected from the discontinuity. 

The forward traveling mode at plane j and the baclward traveling wave at 
plane j+1, 


{A'^h, (33) 

are determined by considering the propagation through uniform sections to the 
left and to the right of the plane of discontinuity. 

The matrix equation relating the solution at plane j-1 to the solution 
at plane j* in the fon^^ard direction, is (referring to Figure 8). 


{a"'-} = [ 






(34) 


where the elements of the uniform section transmission matrix [U] are given 
by 




(33) 


The matrix equation relating the solution at plane j-F2 to the solution at 
plane j-M, in the backward direction, is 




(36) 


where, 


^j-(j+l),-(j+2)^^j^ = 6 


IJ 




(37) 


Similar equations apply to the backward-traveling wave in the left section 
and the forward-traveling wave in the right section. 
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By combining Equations such as (30), (32), (34), and (35), one obtains 
a set of matrix equations which interrelate the solutions at all planes 
(note that each uniform section has two planes, one at each end) between the 
source and duct termination* If, in addition, two equations are written 
establishing the modal source input and reflection matrix at the source 
plane, and the reflection matrix at the termination plane, one obtains a 
complete set of equations for the solution vectors at each plane* By stack- 
ing the solution vectors at each plane into one large overall solution 
vector, a matrix equation is obtained relating the solution vector to the 
source vector in the form 

[S] {A} = {Q} (38) 


where [S] is the system matrix comprised of transmission, reflection, and 
uniform section matrices stacked into the appropriate blocks, {A} is the 
stacked modal solution vector, and {Q} is the stacked modal source vector* 

It is not necessary to locate the source plane at the end of the duct; it 
can be arbitrarily located by stacking the source modal participation into 
the appropriate partition of the source vector* 

The solution of Equation (38) , indicated symbolically by 

{A} = [S]"^ {Q} (39) 

provides the modal participation at each plane, for both forward and backward 
traveling waves. The order of the system of equations to be solved is the 
number of planes multiplied by the number of modes used in the expansion 
(held constant at each plane) times tx^o (for forward and backward traveling 
waves). Complex arithmetic is used throughout. 

Two indicators can be used to gauge the accuracy of the solution. First, 
the acoustic pressure or velocity profiles can be reexpanded at planes 
adjacent to an axial discontinuity to check for continuity. Second, the 
energy flux at adjacent planes can be calculated to check for energy balance. 


2.5 ENERGY FLUX 


The expression for the acoustic energy flux in a duct with uniform sub- 
sonic mean flow (|m|< 1) contains terms representing the no— flow acoustic 
energy flux and convected kinetic and potential acoustic energies. The axial 
component of the acoustic intensity, i.e., the acoustic power flux per unit 
area, is given by (7) 


. U 2 

I =*<pv > + - 2 > 

z ^ z p 

o 


“2 

c 


<pv^> 


P U 


<v®> 

z 


(40) 
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where 


= axial component of acoustic velocity 
U = mean flow velocity in the axial direction 
= density of air 

The brackets <> designate time averages. 


Making use of the momentum equation relating acoustic pressure and velocity 
and considering modal acoustic intensity, one gets expressions for the 
intensity for the cases of 


a. the f onward traveling wave modes 



n) = Re 


p„(p„) 
m n 

2p c 
o 



l+M^+M 



(41) 


b* the backward traveling wave modes 


r(m,n) =Ra 


* 


p^CPf.) 

m n 

2p c 
o 


) + M 

- k k 


(41) 


where m and n are modal indices and the asterisk denotes the complex 
conj ugates . 

The corresponding modal energy fluxes, E^Cm^n) and E (m,n) are obtained 
by integrating the respective modal intensities (41) and (42) over the duct 
cross-sectional area, i.e», 


1 4 

+ /* + 

E^(m,n) =» J I^(m,n)dx 


and 


E (m,n) 
z 


J* l^(m,n)dx 


(43) 


(44) 
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The total energy flux in the forward and backward direction is calculated 
from Equation (43) and (44) by performing double summations over all modes. 

The net energy flux at any cross-sectional duct plane is the algebraic sum 
of the total forward and backward energy flux. 

The analysis detemines the modal energies for forward- and backward- 
traveling waves at the planes of discontinuity, as well as the modal energy 
suras and the net energy flux* This allows investigation of energy transmission 
and reflection effects at duct section interfaces. 


The overall suppression is composed of the combined attenuation of 
individual modes, each at a particular level of participation. The 
individual mode attenuation rates are given by the imaginary part of the 
propagation constant, Equation (5), which can be written as 


K 


m 


a 

m 


-fix 


m 


(45) 


The attenuation rate for each mode is, then, 

Att - -8.686 r H 
m 

where the units are dB per unit duct height* 


(46) 
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SECTION 3.0 


PRELIMINARY OPTIMIZATION STUDIES 


The analysis was applied to optimize the design of a two-element suppres- 
sor for test in a laboratory duct test facility. The method used was to 
deterniine the impedance components of each element, both resistance and 
reactance, which would yield the maximum suppression for the test conditions 
to be imposed in the experiment* 

3.1 SPECIFIED CONDITIONS 


The conditions used for the optimization study were as follows; 

& Rectangular duct geometry as sho\^m schematically in Figure 9 

The height between treated sections was 26.7 cm (10,5 in.); the 
duct width was 10.2 cm (4 in,); the total treated length was 
91*4 cm (36 in.), equally divided for each of Che two treatment 
panels* This resulted in a total treatment length/duct height 
ratio of 3,43. 


& Design frequency of 2000 Hz with input modal energy distribution 
as determined by measurement in the test facility. This included 
both no-flow and flow at 0,3 Mach number. 

This choice of design frequency resulted in a duct height/wave- 
length ratio of 1*535 so that at least two higher-order modes 
were cut-on for the no-flow condition; at least one more was 
cut-on with flovj. Thus, counting the lowest order mode, four modes 
could propagate with flow in the hardwall duct. More could 
propagate depending upon the treatment impedance components 
(to be determined by the analysis). As noted in Section 4*1, 
the sound source was deliberately set off-axis to excite higher 
order modes, 

0 Input modal energy distributions (modal expansion coefficients) 
shown in Figure 10 for no— flow and in Figure 11 for flow at 0.3 
Mach number. These distributions were determined experimentally 
for the untreated, hardwall duct. Note that the distributions 
are substantially different with and without flox>7. Thus, the 
optimized impedance components derived from the analysis include 
more than simple flow convection of the acoustic wave. Also note 
that some energy in cut-off modes is present at the source, but 
the cut-on modes dominate in both cases. 

In the analysis, seven modes were included in the energy balance * 
Continuity of particle displacement was used for the wall boundary* Since 
the tests were conducted with an air supply at room temperature, standard 
atmospheric conditions were assumed. 
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Figure 9, Schematic View of the Rectangular Duct, 
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Modal Indexj m (From Equation d) 


Figure 10. Source Modal Coefficients for 2000 Hz Signal, Mach 0,0, 









3.2 ITERATIVE OPTIMIZATION PROCEDURE 


The analysis was used to estimate the optimum effectiveness first, of a 
uniform liner, and then, of a two— segment liner for the specified geometry, 
frequency, and source conditions. The uniform treatment involves only too 
parameters (the two impedance components), but the two-sectional liner 
requires that four parameters (two impedance components for each liner 
segment) be varied. 

An iterative optimization procedure v/as used, as outlined below: 

0 The best uniform liner was found by parametrically varying its 
impedance components. 

Figures 12 and 13 show the calculated suppression as a function 
of normalized reactance and resistance respectively, for the no- 
flow case; maximum attenuation of about 28.2 dB occurs for Z/pc - 
1.25 - 1.65i.* 

Figure 14 shows contour plots of constant attenuation in the 
uniform liner impedance plane at Mach 0.3 flow; maximum atten- 
uation of about 21,6 dB occurs for Z/pc = 0.9 - l.Oi. This is 
markedly less than that predicted for no-flow conditions. 

0 Using the impedance of the optimized single-element liner for 
the domstream segment, the impedance components (resistance 
and reactance) of the upstream liner segment were varied until 
peak suppression was attained. 

At M = 0, a maximum attenuation of 29.2 dB is predicted (first 
segment impedance ratio Zi/pc = 1.15 - 1.8i, followed by - 
1.25 - l*65i). At flow of Mach 0,3, the maximum attenuation is 
predicted to be 22.4 dB (Zl/pc = 0,9 - 1.2i and Z2/pc - 0.9 - l.Oi), 
The associated contour plots of constant attenuation of the too- 
segraent liner in the impedance plane of the first segment are sho^m 
in Figure 15. 

» The impedance of the upstream liner segment was unchanged while 
the impedance components of the do\mstream segment were varied* 

In this step for the no- flow case, a peak suppression of 32.2 dB 
was calculated (Z^/pc = 1.15 - 1.8i and Z 2 / 0 C = 1.33 - 1.35i), 
Correspondingly, a peak suppression of 25.2 dB was predicted for 
the case with Mach 0*3 flow* (Z]^/pc = 0*9 - l*2i and Z 2 /PC = 0*9 - 
0.81) as shown in Figure 16. 


Note: The normalized impe^an^es given in this report are all referred to the 
more conventional e form. This is the complex conjugate of the 
value which is used in the theoretical analysis. 
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Figure 12, Calculated Suppression lor a Uniform Liner as a Function of 
Liner Reactance Ratio. 
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« The alternating iterative procedure of the preceding too steps was 
repeated until the impedance changes yielded less than one dB 
increase in suppression. 

• Having achieved the latter, the corresponding values of the resis- 
tance and reactance ratios were varied by ± 0.1 each to ascertain 
that the predicted suppression was indeed near or at a maximum. 

3,3 RESULTS QF OPTIMIZATION 

The predicted peak suppressions, and the associated Impedance components, 
at the conclusion of each iteration without and with flow are listed in Tables 
I and II 5 and shown graphically in Figures 17 and 18, respectively. Bar 
graphs of the predicted peak suppressions from each iterative step are sho\ra 
in Figures 19 and 20. 

The final result for the optimized too-segment liner without flow is about 

39.0 dB, an increase of 10*8 dB over the optimized uniform liner suppression 
of 28.2 dB, With flovj the result is 29.3 dB versus 21*6 dB for the uniform 
liner. 

In all cases considered, maximum suppression occurred at negative reac- 
tances of both liner segments. In the course of the iteration for the case 
without flow, the resistance ratio of the first liner segment decreased from 
the single-element optimum value of 1.25 to 0.85, Its reactance also decreased 
from -1*65 to -*1,85 (an increase in absolute value). The resistance and reac- 
tance of the second liner segment increased. The former under^^ent a moderate 
increase from 1.25 to 1.35, while the latter changed more rapidly from -1.65 
to “1.04 (a decrease in absolute value from 1.65 to 1.04). 

For the case with flow of Mach 0.3, the changes in the liner segment 
impedance ratio components followed similar trends as in the case of M = 0 . 

This is shown in Figure 21c 

While . the. resistance, of the first liner segment decreased somewhat (from 
0.9 to 0.7) in the course of the iterations, the absolute value of its reac- 
tance increased (from 1 to 1.37). The resistance of the second segment 
remained unchanged (0.9) , the absolute value of its reactance decreased from 

1.0 to 0*7. 

Although the approximate design procedure of dividing the respective 
resistance and reactance values at M = 0 by the factor (1+M) , with M =0.3, 
does not yield the corresponding exact values of the liner segment resistances 
and reactances at M = 0,3, it produces the expected trends. The 1/ (14M) 
factor would be obtained from the continuity of particle velocity assumption, 
and is also the value obtained from empirical results. Based on continuity 
of particle displacement, a factor of 1 /Cl-fM)^ would be chosen. 

For an optimized two-element treatment arrangement, the final value of 
impedance requires a thinner liner element with a higher' faceplate porosity 
for the first liner segment than for the second one* 
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Table I. Summary of Rectangular Duct Propagation Program Liner Optimization, Mach 0.0, 

2000 Hz M = 0.0 Lj_=L 2 = 45.72cm (18^') L = L^fL 2 H/A. = 1.535 L/H ^ 3,43 

^;l/ P c ^ 2 /^^ 

^ 1^2 H 

^ 

K & Sound 


Iteration 

Number 

Upstream Section 
Impedance Ratio 
Z^/Pc 

Dovmstream Section 

Impedance Ratio 
Z^/Pc 

Total 

AdB 


1 

1.25 - 1.651 

1.25 - l,65i 

-28.2 

(Single Element Optimum) 

2 

1.15 - 1.801 : 

1.25 - 1.651 

-29.2 

(Iterate on 2^/Pc) 

3 

1.15 - i.aoi 

1.33-1.351 

-32.2 

(Iterate on Z^/Pc) 

4 

1.0 - l.S3i 

1.33 - 1.351 

-34.2 

(Iterate on Z^/Pc) 

5 

1.0 - 1,831 

1.35 - 1.151 

-36.2 

(Iterate on Z^/Pc) 

6 

0.93 - 1.851 

1.35 - 1.151 

-37.8 

(Iterate on Z^/Pc) 

6A . 

0.93-1.851 

1.35 - 1.041 

-38.55 

(Iterate both 2j^/Pc and Z^/Pc) 

6B 

0,85 - 1.851 

1.35 - 1.041 

-38.95 

(Iterate on Z„/Pc) 



Table II, Summary of Rectangular Duct Propagation Program Liner Optimization, Mach 0.3. 


£ = 2000 Hz M = 0.3 n 2 L^=L 2 = 45.72ctn (18") L = R/X = 1.535 L/H = 3.43 



M & Sound 


Iteration 

Number 

Upstream Section 
Impedance Ratio 
Zl/Pc 

Downstream Section 
Impedance Ratio 
Z^/P. 

Total 

AdB 


1 

0.9"li 

0.9-11 

-21.6 

(Single Element Optimum) 

2 

0. 9-1.21 

0.9-11 

-22.4 

(Iterate on Z^/pc) 

3 

0.9-1.21 

0,9-0.81 

-25.2 

(Iterate on Z.^lPf^') 

4A 

0.8-1.31 

0.9-0,81 

-26.9 

(Iterate on Z^^/Pc) 

4B 

0,8-1.31 

0.9-0.71 

-28.1 

(Iterate on Z^JP^') 

5A 

0.7-1.31. 

0,8-1.41 

0.9-0.71 

0,9-0.71 

-28.46 

-28,66 

1 (Iterate both Z 2 /PC) 

5B 

0.7-1.371 

0,9-0.71 

-29.31 j 
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Figure 18. Optimum Impedance Components and Associated Suppression in a Li 
Rectanaular Duct with Flow. 
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Figure 22 shows the results of iteration 6B (See Table I). The objective 
of that iteration was to explore the effect of the resistance, particularly 
very low resistance, of the first liner segment on the suppression. Thus, 
the resistance of the first (upstream) segment was varied from very low values 
to high values, while its reactance and the impedance of the second segment 
obtained in the preceding iteration (6 A) were kept constant. 

The second liner segment was found, in all optimized cases, to cause more 
suppression than the first one of equal length. An example at Mach 0.0 is 
shorm in Figure 22. The figure shows, in addition to the total suppression (due 
to both liner elements), the attenuation due to the first segment only. Tiie 
latter attenuation is close to zero for near-zero resistances of the first 
liner segment; it increases with resistance and attains a broad maximum for a 
resistance ratio of 0.85. Beyond that, it decreases very slowly. The total 
attenuation of both segments peaks at the same resistance as the first segment. 
Its increase and decrease, however, are characterized by steeper slopes than 
the attenuation by the first liner segment. 

/ Q \ 

For very similar conditions , Baumeister^ has found optimized two-element 
configurations for very low resistances of the first section. During the con- 
vergence process in the search for the optimum, special effort was made to 
check for maximum suppression conditions which decreased the first liner 
resistance. There was no indication of any trend toward low first liner 
resistances In this case. In fact, all indications were that the optimum 
suppression occurs at a unique value of impedance (holding the panel 
lengths, and the modal content constant), the particular optimum cases 

run by Baumeister were attempted by the modal analysis, numerical convergence 
problems were encountered in most cases, leading to indeterminate results. 

As will be shown in a later section, there are optimization conditions (in 
particular, higher n values) which lead to low resistances in the first 
section. 
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SECTION 4,0 


DUCT TEST PROGRAM 


The test program had as its primary objective the experimental demonstra- 
tion of the suppression properties of the optimized treatment designs deter- 
mined analytically in Task I of the program. The test program was accomplished 
in two phases. Phase I comprised the measurement of suppression for six 
treatment configurations using liner design parameters which were believed to 
bracket the desired impedance compouents. Phase II ti/as a test of several 
treatment configurations with the objective of improving on the results of 
Phase using knowledge gained from the first series of tests. 


4,1 DUCT FACILITY PESGRIPTIQN 

The rectangular cold flow duct facility provides the capability of measur- 
ing treatment sample acoustic transmission losses under conditions of continu- 
ous flow with Mach numbers up to about 0,5. The duct is 10*2 cm (4 inches) in 
width with har dwell sides and allows duct heights up to 40,6 cm (16 inches) 
with treatment panels at the top and bottom. The treatment panels are con- 
tained in trays 91,4 cm (36 inches) in length which allows interchangeability 
of perforated faceplate and honeycomb backing. The treatment panels can be 
segmented in the trays, allowing axial variation of faceplate porosity and 
cavity depth in the treatment section. 

The noise source for the duct is a Ling Model LPTlOO electropneumatie 
high intensity driver which can be excited by broadband or pure-tone input 
signals. As shown in Figure 23, the source is mounted upstream of the treat- 
ment, providing an exhaust diiet mode of operation. Note that the source input 
port to the duet is mounted asymmetrically with respect to the duct vertical 
centerline for thase tests* This was done specifically to insure the presence 
of higher order duet modes in the treatment input signal. 

The duet terminates downstream of the treatment with an unflanged. end, 
exhausting into a reverberant room, A muffler section is included at the end 
of the duct to minimize end reflections and waves reentering the end of the 
duct. 

The acoustic measurements are taken with traversing probes iJistalied up- 
stream and downstream of the. treated section* Tte probes traverse vertically 
across the duct along a line midway between the sidewalls of the duct. The 
probes are specially construeted using a dynamic pressure transducer that is 
a 0,23 cm Kullte sensor mounted in the tip of the probe, as shown in Figure 
24, The Kullte transducer is a highly linear device based on a straiti gauge 
principle* It has lower sensitivity than the more conventional B&K condenser 
microphone but iiS much less sensitive to meGhanicai vibratioh* The low level 
output of the Kiilite transducer requires the use of a Ibw-^noise^ high-gain 
preamplifier for signal conditio niiig* Princeton Applies Research Model 113 
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Figure 23. Riectangular Duct Test Configuration. 









preamplifiers were used for eacb. Kulite. The tip^mouuted Kulite probes were 
used for both continuous SPL traverses and modal nieasurements. 

The modal measurements required the determination of the cross-spectral 
density betv 7 een the probe and a reference micronhone as a function of probe 
immersion across the duct (see Section 2.3), The reference microphones were 
B&K 1/4-inch (0*635 cm) microphones mounted flush in the sidewall of the duct 
at the same axial positions as the traversing probes (see Figure 23)* The 
cross“spectrum was found using a Princeton Applied Research Model lOlA Corre- 
lator and Model 102 Fourier Analyzer. These instruments can be set to read 
out the magnitude and phase of Che cross-spectrum at a single frequency while 
the probe is traversed across the duct* Using a 2n-second calculation time 
constant in the analyzer, the probe was slowly traversed across the duct as 
the analyzer continuously updated the calculation* This method allows a 
continuous plot of the cross-sectrum (and thus complex acoustic pressure) 
profile to be obtained at the frequency of interest. 


4*2 PRELIMINARY 


The analytical optimization study requires the modal participation at the 
designated source plane as input to the computation. Before initiation of the 
analytical study, a modal measurement was performed in the 26.7 cm high duct 
in a hardwall configuration to supply this input. It was assumed Ghat the 
source characteristics would not change appreciably for the treatment config- 
urations to be tested later* 

Figure 25 shows the complex pressure profile obtained at 2000 Pz and 
Mach 0*3 in the duct. Figure 26 presents the modal decomposition of this 
signal in hardwall rectangular duct modes* Figures 27 and 28 show the com- 
plex profile and modal expansion for the signal at 2000 Hz and Mach 0*0, which 
was also measured at this time* 


4,3 PHASE I TE8T PROORA}! 

4*3.1 Obj ectives 

The obiective of this first series of tests was to demonstrate the opti- 
mized two-element treatment design found by analysis. The analytical results 
indicated that attenuation would be maximized if the two treatment sections 
were given the following impedances (see footnote in Section 3.2): 

Section 1 (Upstream) Z/pc =0*7-1.41 

Section 2 (Do™stream) Z/pc = O.q - 0.71 

This design is based on a source characterized bv the modal participation 
given by Figure 26, at a mean flow of Mach 0.3 and a frequency of 2000 Hz. 

The predicted attenuation in PT-JL under these conditions is 29.3 dB* 
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Figure 25. Complex Pressure Proltle, Upstream Vvnhe. ^ Prelim- 
inary Experiment, Mach 0,3, 2000 Hz. 




Figure 26. Modal Participation. Preliminary Experi- 
ment, Mach 0.3, 2000 Hz. (Repeat of Figure 11)* 
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; 27. Complex Pressure Profile, Upstream Probe, Preliminary 
Experiment, Mach 0,0^ 2000 Hz. 
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Figure 28. Modal Participation, Preliminary Experi- 
ment, Mach 0.0, 2000 Hz. (Repeat of Figure 10), 
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4*3.2 Treatment Design and Test Program 


The type of treatment to achieve the optimum impedance components was 
chosen to be single-degree-of-f reedom (SDOi") perforated plate ~ honeycomb 
cell sandwich, panel . The acoustic properties of these SDOF panels have been 
established^ * "^design parameters (faceplate porosity, thickness, hole 
diameter, honeycomb cavity depth) to the wail impedance. Preliminary design 
of the treatment was made using the following expression for resistance 


R/pc 


Q.3M 

a 


(47) 


where M is duct Mach number and c is faceplate porosity. The following 
expression was used for the reactance 


X/pc 


ftp + 


cot ki 


(48) 


where tp is faceplate thickness, u> is circular frequency, k is the wavenumber, 
i is the cavity depth, and 6 is an empirical end correction depending upon 
hole diameter and Mach number* Using results of these impedance calcu- 

lations as a guide, the six configurations listed in Table ITT were chosen 
to he tested . 

For each of the configurations, APWL measurements were made for third- 
octave bandwidths from 500 to 10,000 Hz and for a 50 Hz bandwidth narrowband at 
2000 Hz. For the third-octave measurements, the duct was excited by a high 
intensity source signal which was filtered in a third-octave bandwidth at 800 
Hz* In this case, the higher frequency noise t»ccurs due to rich higher harmonic 
content of the 800 Hz signal. For the 2000 Hz narrowband, the input signal was 
broadband noise filtered in a third-octave bandwidth about 2000 Hz with a super- 
improved 2000 Hz pure tone. The treatment configuration demonstrating the 
greatest APWL for the 2000 Hz narrowband was subjected to detailed modal 
measurements for correlation with analysis. 


4.3*3 Test Results 


Figures 29 to 34 are plots of the measured tl. i rd-uctave suppres- 

sions for the six Phase I test conf igu rat ions , respectively, ior Mach 0.0, 
0.2, 0.3, and 0.4* Table IV lists the 2000 Hz narrowband IHUl. for tiie 
six configurations* 

Examination of the measured third-octave suppression curves in Figures 
29 to 34 indicates that in most cases a higher suppression was obtained at 
2500 Hz than at 2000 Hz* TFiis can be explained by the presence of a strongly 
excited higher order mode whJxh. becomes cut-on in the center of the 2500 Hz 
band* The treatment is very effective in attenuating this particular mode 
even though it is designed for 2000 Hz, 
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Table III* Phase I Test Configuration Definition* 


SECTION 1 SECTION 2 

Porosity Cavity Depth Porosity Cavity Depth 


Config. 0 

HARDWALL 





HARDWALL 




Config, 1 

147. 

1.5 

cm 

<0.6 

In.) 

107. 

2.3 

cm 

(0.9 

In.) 

Config. 2 

147. 

1.5 

cm 

(0.6 

In.) 

7.57. 

1.5 

cm 

(0.6 

In.) 

Config. 3 

107. 

1.3 

cm 

(0.5 

In.) 

7.57. 

1.5 

cm 

(0.6 

In.) 

Config. 4 

107. 

1.3 

cm 

(0.5 

In. ) 

107. 

2.3 

cm 

(0.9 

In.) 

Config. 5 

147. 

1.5 

cm 

(0.6 

In.) 

147. 

2.3 

cm 

(0.9 

In.) 

Config. 6 

22.77. 

1.5 

cm 

(0.6 

In.) 

107. 

2.3 

em 

(0.9 

In.) 


All iaceplate thicknesses „08 cm, all hole diameters .16 cm 



Corrected Transmission Loss, dB 





Corrected x^ansoiission Loss, dB 
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Corrected 'iranstnis 



Figvire 31. Measured Third-Octave Transmission Loss, Con iigurat ion 3. 
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Corrected Ij^ansmls 



Frequency, kHz 


Figure 32* Measured Third--Octave Transmission Loss, Configuration 4, 
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Corrected ‘iransmission Loss, dB 



Figure 33* Measured Thlrd-^Octave Transmission Loss, Configuration 5, 
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Frequency ^ kHz 


Figure 34* Measured Third--Octave Transmission Loss, Configuration 6 






Table IV. Phase I Test Results -* Measured Narrowband Transmission Loss 
at Macii 0.3, 2000 Hz, 




SECTION 1 

SECTION 2 

AlWL, dB 

Conflg. 

1 

(147o> 1*5 em) 

(10%, 2.3 cm) 

-17,5 

Configo 

2 

(14%, 1.5 cm) 

(7.5%, 1,5 cm) 

-11.0 

Gonfig. 

3 

(10%, 1,3 cm) 

(7.5%, 1.5 cm) 

-12,5 

Conflg* 

4 

(10%, 1.3 cm) 

(10%, 2,3 cm) 

-19.0 

Conflg, 

5 

(14%, 1,5 cm) 

(14%, 2.3 Cm) 

-17*0 

Conflg o 

6 

(22,7%, 1,5 cm) 

(10%, 2,3 cm) 

-13.5 

Conflg, 

(Repeat) 

4 

(10%, 1.3 cm) 

(10%, 2.3 cm) 

-21.5 


In Figure 35 an empirically determined rectangular duct suppression 

prediction curve j based bn tlie results of a large number of 

single-element liners, is compared to the measured third-octave suppression 
for Gpnfiguratldn 4* Also plotted bn the Mach 0.3 graph are the 2000 Hz 
narrowband and third-octave measurements made with the 2000 Hz third-octave 
source , which, as noted above, was different from the source used for the 
rest of the third-octave measurements. Note that the different source 
characteristics at 2000 Hz have given a difference in measured suppression 
of 3*5 dB* the narrowband measurement,, made at only one frequency, shows an in 
crease of 5*5 dB over the third-octace with the same source characteristics* 

At Mach 0.0 and 0 . 2 , the measured sappressiqn agrees quite clpselv 
x^ith the empirical prediction curves . for single-element liners . Ah the 
design Mach number of 0 . 3 , the tx7o “element measured suppression using the 
2000 Hz third-Qctave source characteristlG (more representative of design 
GonditiPiis) exceeds the empirical predictib^^ for the best single phabe 
liner by about 3>5 dB . For third'^pctave- frequencies above 2000 Hz, the 
measured suppression exceeds the empirical prediction by: a substantial margin. 
At Mach 0.4, the high frequency enhancement of the two- clement liner is still 
in evidence, although ho improvement over the empirical data is noted at 
2000 Hz. 

One can conclude that in this case the two-element liner docs not appear 
to offer a significant advantage over the single-element liners for frequencies 
lower than the design . frequencyV but provides signif lean t.Ty enhanced suppress - 
Sibil for higher frequencies at Mach numbers equal to or greater than the 
design Mach number. Note that the narrowband suppression peak within the 
2000 Hz third-octave, at the design condition, is very sharp x^ith respect to 
the overall third -octave band suppress ibn% 

Con f igurat ion 4, wh ich came elo s est to the predicted op timum suppre s s ion 
of 29.3 dB, x^as chosen for detailed modal measurements . A repeat narrox^band 
APWL measurement of this ^ configuration a suppression, of 21.5 dS. overall, 

(compared with 19*0 dB previously) and 8 dB over the first section of treats 
ment only. The complex pressure profile at the upstream probe Ipcat ion is 
shown in Figure 36, and its expansion into duct modes in Figure 37. This 
; variation' -in mbabureTheiit Wab • apparbnily •' catised-; -"by •''•slight" dif fefen&es ' "in ^mSdal • 
content, as evidenced by somexi^hat different SPL traverse profiles for the 
two cases. In subsequent testing it was found that slight differences iti 
modal content can be expected in the duct with each different test case, and 
that these small differences in imbcial bbrltent dpuld produce AP 
of this magnitude. This effect is considered in more detail below. 



Corrected ikansmlsslon Loss, dB 


Measured 50 Hz Narrowband Transmission Loss 
Measured Third- Octave Transmissioa Loss Using 
2000 Hz Source 

Empirically Determined Third-Octave Suppression Predi 
Measured Third-Octave Transmission Loss Using 
SOO Hz Source 


M=0.2 ■ 




Frequency^ kHz 


Figure 35* Measured Third^Octave Transmission Loss^ Configuration 4 
Compared to Empirical Prediction for Rectangular Duct. 

Measured Thlrd*-Octave Data from Figure 32, 












4.3*4 Analytical Correlations 


Based on the modal measurements taken at the plane of the upstream probe, 
an analytical prediction of the suppression was made* For these calculations 
at Mach 0*3, the liner impedances were taken to be (Equations 45 and 46): 

Section 1 (10% porosity, 1,3 cm deep), = 0.9 - 1.4i 

Section 2 (10% porosity, 2,3 cm deep), ^2^^^ = 0.9 - 0.71 

The results are: 

Measured Calculated 

APWL in Section 1 8 dB 6.0 dB 

Overall 19 to 20.6 dB 

21.5 dB 

These results indicate that the suppression was predicted quite closely 
for the given modal input and liner impedance* Neither the modal input nor 
the liner impedance is exactly the same as the conditions used for the opti“ 
mization analysis, which explains why the predicted 29.3 dB for the optimized 
case was not reached. Figure 38 shows a comparison of the modal content for 

the original set of optimization modes and those in the input signal to con- 

figuration 4, This variation in modal content is discussed below* Note for 
reference that the predicted uniform liner optimum suppression was 21*6 dB. 


4.4 PHASE II DUCT TESTS 


4*4,1 Phase II Test Objectives and ReoptlmlzatlQn of Design 

The results of the Phase I testing indicated tiiat the optimization condi- 
tions were not met in respects; first* the modal content had changed, 
and second, the impedances differed from the optimum values* The objective 
for the second phase of the test program was to allow a second Iteration at 
meeting the optimum design conditions. It was decided to reoptimize the 
design of the two“element liner based on the new modal content measured in 
Configuration 4 of the Phase I tests. 

Using the modal input defined in Figure 37, the following design was 
found to maximize suppression for the two-segment liner with equal length 
segments at 2000 Hz, tlach 0.3, 

Section 1 Impedance, = 0*6 - 1.2i 

Section 2 Impedance, Z^! Qc - 0*7 - 0.5i 

Overall Attenuation APl^L = -35*9 dB 
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Figure 38. Comparison ol Measured Modal Participation, Preliminary Experi- 
ment and Configuration 4, Mach 0.3. 2000 Hz. 
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Note that the impedance values are only slightly changed from the previous 
optimization case, but the predicted suppression has increased by over 7 dB- 
This variation is due to the difference in modal content between the two cases. 

In addition to testing a treatment configuration to attempt to meet the 
revised optimization conditions, it was decided that other objectives of Phrse 
II testing would be to investigate the repeatibility of the modal pattern and 
the sensitivity of the predicted attenuation to uncertainty of the impedance 
components. If the modal pattern varies appreciably with the treatment con- 
figuration, or is highly sensitive to frequency, it would be extremely diffi- 
cult to obtain the optimum condition, particularly in view of the sensitivity 
of the peak suppression to the impedance components. The sensitivity to 
impedance changes about the optimum value is indicated in Table V in which 
APWL*s have been calculated for cases where each of the four specific impe- 
dance components has been varied independently by ±0.1 from the optimum value. 
Using the analytical models for determining the impedance components introduces 
an unkno™ amount of error relative to the actual impedance of the liners in 
the duct. Tlie sensitivity of the predicted suppression to expected errors in 
impedance component determination was examined. 


4,4,2 Phase II Test Configurations 

For the Phase II duct tests, the two-element liner was designed to attempt 
to meet the revised optimization criteria. Based on the analytical impedance 
models, the following liner parameters were chosen: 

Section 1 14% porosity, 1.8 cm (0,7 inch) cavity depth 

Section 2 12% porosity, 2,5 cm (i.O inch) cavity depth 

Both sections of this configuration, denoted Configuration 7, had 0.08 cm 
(0.032 inch) faceplate thickness and 0.16 cm (0.0625 inch) hole diameter. In 
addition to this configuration, three other configurations were run as part 
of Phase II testing. The other configurations consisted of Configuration 7 
with the order of the segments reversed, and two cases with the same cavity 
depths as Configuration 7, but with higher and lower porosity of the face- 
plates , respectively . 


4.4.3 Phase II Test Resul ts 

A list of the Phase II test configurations with the measured narrowband 
suppressions for each case is presented in Table VI. The results of the 
source modal measurements for Configurations 7 to 10 are present in Figures 
39 to 46. In evaluating the test results which follow, it is important to 
note the variation in modal content which occurred from one test configuration 
to the next. This variation was more than initially hoped, and made achieving 
the optimum conditions difficult. 

Before measuring transmission loss at 2000 Hz for Configuration 7, 
it was decided to investigate the sensitivity of the tuning to frequency in a 
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Table V. Calculated Suppression for Two^Element Liner, 



Configuration 4. 



^ 2000 Hz 

• Mach 0.3 


IMPEDANCE 

IMPEDANCE 

aPWL ♦ dB 

Section 1 

Section 2 


«6 -1.2i 

.7 -c5i 

-35.9 

«5 -l»2i 

-7 “„5i 

-30.1 

o7 -1.2i 

„7 -.5i 

-27.5 

.6 -l.li 

,7 -„5i 

-29.9 

*6 -1.3i 

«7 -.5i 

-27.5 

*6 -1.2i 

„8 -„5i 

-31.9 

b 6 -lo2i 

u6 -«5i 

-31.4 

«6 -lu2i 

-7 -.4i 

-30.1 

„6 -1.2i 

„7 -«6i 

-34.7 


NOTE: Modal input given by Figure 37. 



Tfitjlr VI, 11 Tr;, t. Con f'i t i on Dei 3 nation 


J h'C i 1 oil 

Con 1 1 t i on 

1*0 1 1 t y 

% 

Cnv i l y 

Dep t h 

Nominal 
Inipeciance 
Z/q c 

cm 

1 n. 


7 

1 1 

! .H 

0.7 

0.64 - 1.15i 


H 


2.5 

l.o 

0.75 - 0.55i 

\ 

H 

:^o 

1 .s 

0.7 

0.43 - 1.301 




1 .H 

0.7 

0.75 - 1.201 


7 

12 

2 . 5 

1.0 

0.75 - 0.551 


.s 

1 1 

l.H 

0.7 

0,64 - 1.151 


n 

H 

2,5 

1 .0 

0.64 - 0.601 


m 

lo 

2.5 

1 .0 

0.90 - 0.551 


|u^ I ;> jx\*\ I t!' Kmv'i •]> t ;i 1 r s i T^(' 0,US cm Thick* CK 16 cm Hole 
n \ ;itnc 1 1 - t 


;mv: imu’ciM.ur 0,0S cm Thick, 0,08 cm Hole Dlnmeter 



















range about 2000 Hz, To do this, an automatic level control device was con'- 
nected to the input of the Ling driver such that the sound level at the input 
plane probe to the duct could be maintained at a constant level as input 
frequency was slowly swept from 1450 to 2350 Hz, using the upstream probe 
signal as a feedhatk control. With both probes in the fully retracted posi- 
tion (near the lower wall) , the SPL signals from both probes were recorded in 
20 Hz bandwidth during the frequency sweep. 


Although the actual APWL cannot be inferred from ASPL's at one immersion, 
from previous data it was felt that the ASPL at the lower wall would be repre- 
sentative of the APWL, thus giving an indication of the best tuning frequency* 
Figure 47 is the result of this frequency sweep. Note that a high ASPL 
occurs at about 1950 Hz, but the ASPL right at 2000 Hz is much less. From 
this, it was decided to take the AP\^ measurements in e narrowband about 
1950 Hz for this configuration. In this and the following cases, the source 
was excited with a pure- tone and the measurement is in 20 Hz bandwidths. 

Figure 48 is the measured narrowband SPL traverse for Configuration 7 
for the upstream and downstream probes. The APWL based on these measurements 
is “^20 db. This falls short of the predicted optimum attenuation by about 
16 dB. 

Using the modal input given for Configuration 7 by Figure 40, and the 
analytically determined imped ances of 

H = 0,64 - l.lSi 

PG 

fg = 0.75 - 0.55i 

pc 


the suppression predicted for Configuration 7 is 28.3 dB- This overpredicts 
the measured suppression by over 8 dB. The most likely cause of this discre- 
pancy between the measured and predicted suppression is the error in deter- 
mining the impedance. An indicatioii of the sensitivity of the analytical 
program to variations in impedance about the optimum value has already been 
indicated in Table V. The nominal value of the impedances of the liners 
in Configuration 7 are all within O.OSpc of this optimum value, and are pre- 
dicted to lose only 7.6 dB in suppression from the optimum value, even with 
the different modal content of Configuration 7, compared to that used for 
the optimum. 

Since there is no way of determining precisely what the liupcdance com*^ 
ponents in the duct actually are, a study was made of the expected variation 
in the impedance components as functions of the independent variables in the 
analytical impedance expressions, and of the sensitivity of predicted suppres- 
sion to typical variations o^i the impedance parameters. This error analysis 
also gives an indication of the practicality of designing a treatment panel 
to be optimized for the suppressiGn of a ptire-tone or narrowband signal. 
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Figure ^ Ne^rowbari^i; (20' SPL Traverses, Configuration 7, Macli 0. 3, 1950 Hz j ApWL =: 20 dB, 



By Equation 45^ the signif iGant variables for the resistance model 
are the Mach number M and the porosity a. Figures 49 and 50 show the 
variation of R/pc for the two liners over the expected variation in M and a. 

The reactance depends upon the porosity; a, the cavity depth; £, the faceplate 
thickness; tp , the frequency; and the end correction 6, which in turn depends 
on the faceplate hole diameter and Mach number* Several of these variables 
can be eliminated as having relatively small effect. Figure 51 shows the 
empirical curve of the coefficient used in the end correction term as a func- 
tion of M* ^Tote that at Mach 0.3 6 is nearly zero and can be ignored. The 
error due to uncertainties in faceplate thickness should likewise be small. 

The variations of reactance for the 1*8 cm (0.7 inch) and 2.5 cm (1.0 inch) 
cavity depth liners due to expected variations in M, a, and to are shown in 
Figures 52 and 53, respectively. Combining the largest of these variations 
for the two liners leads to rough ’^expected error” plots shown in Figures 
54 and 55* Based on these plots, the expected variation in impedance for 
each impedance component for the suppression sensitivity study is taken to 
be tO.Zpc. 

The variation in predicted suppression for Configuration 7 as a function 
of the variation in impedance components about the nominal values of 

= 0.64 - 1.15i 
pc 

^ = 0.75 - 0.55i 
pc 

is presented in Figures 56 — 39 and in Table VII. The figures show the 

variation as the impedance components are varied independently, holding the 
other three constant. The table shows the results of varying a number of the 
components in various combinations, as noted* 

The results show the suppression in this case to be particularly sensi- 
tive to the resistance of the first liner and the reactc^nces of both liners. 

The values of predicted suppression obtained vary from a low of 15*7 dB for 
the case where both resistances are lov; and both reactances less negative 
to a high of 35,4 dB for the case where X^/pc is +D.2 from the nominal value 
and the rest of the components are nominal* Note that certain sets of imped- 
ances can be found i^hich give higher suppressions than the "optimum” impedance 
for this case, which vJas based on the different modal content of Configuration 
4, The most important ccnelusion from these results is that the measured 
suppression of 20 dB is within the range of expected error in the impedance 
comporLents, and thi.s is the moat likely cause of disagreement in the predicted 
and measured values* 


For Configuration: 8, which is Configuration 7 with the panels reversed 
in axial position,, the impedances used for calculation at 1950 Hz were 

Zj^/pc = 0.75 ■* 0.551 Z^/pc = 0.64 - 1.15i 
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Figure 53. 


Variation of Reactance Ratio as a Funet 
Frequency for a 2.5 cm in.) Nomina 













Figure 54, Expected Error in Determination oi Impedance 

Components j 14% Porosity, 1.8 em (0,7 in*) Deep 
Liner* 



Figure 55* Expected Error in Determination of Impedance 

Components, 12% Porosity, 2*5 cm (1*0 in.) Deep 
Liner* 
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Figure 56, 


Variation of ZiPWI. as a 
Function oi R^/Pc r 
Coniigu ration 7, 1950 Hz. 


Figure 57, Variation of APWL as a 
Fuinetion of Xi/pG, 
Configuration 7, 1950 Hz 











Table VII. Variation of Predicted APWL for Configuration 7 with Combinations of Impedance Goraponents 


DESCRIPTION OF VARIATION 

IMPEDANCE OF 
FIRST SECTION 

Zj^/Pc 

IMPEDANCE OF 
SECOND SECTION 
Z^lPc 

PREDICTED 

AFUra. 

dB 

Nominal Case 

.64 - 1 . 15 i 

.75 

-. 55 i 

- 28.3 

Both resiS'Zances low 

.40 - 1.15 

.55 

-.55 

- 17.9 

Bo til resistances high 

.80 - 1.15 

.95 

-.55 

- 26.2 

Less negative reactances 

.64 - 0.95 

.75 

-.35 

- 25.4 

More negative reactances 

.64 - 1.35 

.75 

-.75 

- 19.7 

Both resistances loW;, both 
reactances less negative 

.40 - 0.95 

.55 

-.35 

- 17.0 

Both resistances high, both 
reactances more negative 

.80 - 1.35 

.95 

-.75 

- 19.5 

Both resistances low, both 
reactances more negative 

.40 - 1.35 

.55 

- .75 

- 15.7 

Both resistances high, both 
reactances less negative 

.80 - 0.95 

.95 

-.35 

- 25.2 



Using the measured modes from Figure 42 gave a predicted suppression of 
17.7 dB, compared with a measured value (with the nominal impedances listed 
in Table VI) of 11.0 dB, Using a slightly higher restance of 0.4M/a for 
comparison gives a predicted suppression of 15,4 dB in this case, overpre-* 
dieting by 4,4 dB, Thus, reversing the order of the liners caused a 
substantial reduction in suppression. 

The suppression for Configuration 9, with higher porosity liners, was 
measured at 1940 Hz. For the prediction in this case, the impedance components 
were taken to be 

Z^l pc = 0.43 - 1.3i 
Z 2 /PG = 0.64 ^ 0.6i 

Using the modes from Figure 44, the suppression is predicted to be 
17.6 dB, which compares favorably with the measured value of 16.5 dB. 

For Configuration 10, the porosities of the liners were decreased. In 
this case, to better define the variation in suppression around 2000 Hz, a 
series of narrowband APWL measurements were made every 20 Rz from 1880 Hz to 
2100 Hz, using pure tone excitation and 20 Hz bandwidth data reduction. 

The results of these measurements are shown in Figure 60. The 26.2 dB 
suppression measured at 1900 Hz is the highest measured so far, and was chosen 
for modal measurement, as shown in Figure 46. In addition, a repeat of the 
narrowband APWL measurement gave a value of 26.1 dB suppression. 

The nominal value of the impedance in this case is 

Z^/pc = 0.75 - 1.2i 

Z 2 /PC = 0.9 - 0.551 

This results in a predicted attenuation of 33.5 dB, overpredicting the mea- 
sured value of 26-2 dB by 7.3 dB. It is suspected that the same problems with 
errors in impedance occur here as for Configuration 7. 

Third-octave transmission loss measurements were also made for Configuration 
IG. The results are shown in Figure 61, compared to the empirical duct test 
prediction for single-element liners. This configuration shows a slight im- 
provement in third-octave suppression over Configuration 4 (Figure 32) at 
2000 Hz. The enhancement over the single-element liner prediction, however, 
occurs mainly for frequencies below 2000 Hz. This is most likely to be caused 
by the slight Increase in cavity depths for Configuration 10 compared to 
Configuration 4, which would tend to lower the third-ectave tuning frequency. 
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4.4.4 Modal Explattatlpn of pifffereiice la Suppression 

Further insight to the process of modal attenuation can be gained by 
examining the details of the modal solutions £ 01 : the e^qjerimental Gorrelatipn 
cases. Configurations 4, 7, 8, 9, and 10. Figure 62 shown the plane numbering 
system used ^'or the analytical model. For convenience, the upstreaij treatment 
section will be referred tb as the "first sectibn'*, and the downstream 
treated section as the "second seetlon". 

Figures 63 through 72 present tahulations of the ma^itude and phases of 
the holutton ifibdal GGefficients for the fbrx<rard and bacldjard waves at. each 
plane, the order of the coefficients can be identified with the following 
modes at each plane: 

1 - First symmetric 

2 - First antisyrametric 

3 - Second symmetric 

4 - Second antisymmetric 

5 - Third symmetric 

6 - Third antisymmetric 

7 - Forth symmetric 

Plotted alongside the modal coefficients at each plane are the (analytical) 
complex acoustic pressure profiles at each plane. 

Table VIII is a list of the attenuations of individual modes in the 
first and second sections of each configuration. Only the first four modes 
are included . Table IX is a list of the forward and backward traveling 
energy at each plane and the suppression of forward traveling energy in the 
first and second sections for each configuration. . 

The differences in suppression among each of these configurations can be 
explained by considering: 

a. The relative modal weighting at the source plane, 

b> The modal redistributions in the forward traveliag wave between 
Flanes 2 and 3 and Planes 5 and 6,. 

c. The Individual modal attenuation rates in Sections 1 and 2, v.v - 

d. ' The bverhll energy attenuations in Sections .1 and 2. 

\>/hile some of the effects have obvious explanations, others are the result of 
a very subtle combination of modal redistr-ibutibn and: attenuatiQn. properties . 

At the source of Configuration 4, most bf the energy is in modes Ij 2, 
and 4, with 1 dominating. Passltig into the first sectioh., at plane 3, the 
ennrgy ■ is redistributed predominatly into : mode , S * -Mode 3 is - st.i:W dpniitiaTit 
at the eri^ of the first section', and f einains: so at the end of ' the Second 
section at plane 6. -'f;:;- . '/ 
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Figure 63. Modal Participation in First Seven Modes at Each Plane, 
Configuration 4. 
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Figure 65, Modal Participation in First Seven Modes at Each Plane, 
Configuration 7, 
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7.94i,33bt-ll 

b,ab9oft3£ Cl 

2.70573tP-a3 

«I M^lTuBE oa 

0. 

dr 





3.12a939F-02 

rr5*ab098«F 01 

«^T7B37lb-dl ' 

l,375a39F oa 

t.97l7Biir-g2 

l,9070h^b 01 

fi*47^2bbE.-07 

ra 


i.mflaoiF 02 

I^9ia3?bt-U 

-1 ,aob{j8'jt. oa 

ti 5b?a9ttP-02 

t,94ft37b£ ai 

a^7'j06l^b“Di 

l*uib2abP fra 

S,T796l9fc-p7“ 

-b,70ia2'i£ Q\ 

a*j77h23b-Dl 

-b,3aai7f<p 01 

5.5a91«a£-|l 

&,Bi9b57LJJ0 

3.8fafad82F“02 

-a.C«9i>3Bt 01 

EJtt27a7£-U 

0* 

5^399033fc-03 

-U7bub3UE 02 

g^75lBfaqg-D2 

-1 pfl07 177E, oa 

1 .232591E-Q2 

-i.oio9?bt oa 

3. 6 59 662 F- 02 

-l,ai767lL 01 

q*724l7yNQ3 

6,511Q43£ 01 

S,^0?Jj5frt-92 

-.J4,B5237?E Q{ 

a.«OBttl9F-j3„ 

0, 

S*7«ea9lt-02 

-‘l,71bi7faE 02 

0. 

di 

9*91021frfc«03 

1 .570398E 02 ._ 

3.aU50fr3£-06_ 

I .7S74g2£ oa 

a",0332l!>E-05 

-l,7b‘?^f,bE oa 

1 *149h65E-10 

-l.baieBbb 02 

2,flR2957t-0^ 

-,3,98'3ti?9h 00 

0, 

0. . 


l*2aa95iE-02 

1^69o/77E 02 

I,44fa533t-fl2 

B,b225bib 01 

4.39oe4at-oi 

r,3,ib097B£ 01 

!^l97b?4F-02 .. . . 

1.0394B4E 02 

2*6b023bt-a3 

-a,l497R2E 00 

4 59a27eF-09 

01 

i,2>ai kb r_- Q 3 

=d.dH33a^b_02 

?.l2?7tPF-0i 

7^fe4l3?6b 01 

1.67717SF-0b 

•«a.tJj08bb 00 

■ ' 3*a5tMJ4bt-03 

-a,B97393b 01 

9.a5020lF-ll 

01 ._ 

9^4bflb7bF-04 

-b.SlilbtrE 01 .. 

7.275955E«tl 

*l,2bBb99b 02 

a,0l9632E-03 

l^0U5394t 02 

t,37oa«6E*^a2 

-l,7bO«33E oa 

0. 

tu 

2,397b£9£-03 

i,babb7ab oa 

0* 

Ot 

l*'5bab7bE-02 

-4,749*#72t 01 

0, 


9*5io7Pht-rj3 

-.b,3b7b2bE Ot 

0* 

0, 

2*?035?7f.-0i 

-i,ovsub£ oa 

0 ^ 

Ot 

4,^l«b37Me-04 

tt.OlOFjttb 00 

^0* 

0, 

a.4ba9afaE-ou 

»9 ,3BHtlij£iE 01 

O' .. 

^ 


l,370a«6E-n2 

t,,05B427£ 01 

0, 

0. 

2r397bb9t-fl3 

a,lDbSH3?t 01 .. 

Q j 

0, 

l.5bPh76r-aa 

I,u99bnlb 02 

0^ 

0, 

9*5707B<tr-Oi 

->i.l4?7>i7L 01 

0 J 

_iu 

a,69b66Db-Ob 

l^tffbUPCt 02 

0. 

di 

2^8b7au9r-07 

-9th972bub 01 

0 

- dg. 

b,492179b-Qtt 

L,6bl2B3£ 02 

0, 

0. 


69, Modal Participation in First Seven Modes at Each Plane, 
Configuration 9, 


..-■T- -rl L- .'"a Hr ... 

r", - ■ 


95 















PLANE 


I 


2 


3 


4 


5 


6 


7 


8 


Figure 


FORWARiMfe^E- — -BAC-Kt-7:^SB-l-?AVE 


HiCNITUDL 

PHAS£ 

hacnitude 

PHASF 

2;a60ia7E 00 

-9,2til5fur 01 

5;S0i*S85E-aa 

1,6373o6E 02 

00 

isXUmiE 01 

4,037fa?9F-Ol 

-S.qaoBbOE 00 

i;3S7536E 00 

t,tf22ntJF 02 

!,52c7S3F-01 

1,5553oO£ 00 

2>2707<.E 00 

^l,20?7lo£: 02 

00 

I,136BlBE 02 

7,3OO7?0F-OI 

a*920qb3E 01 

a.obo44ot^o3 

1,5B5294E 02 

5*S507?flF-Ol 


5*3l73fi6E-03 _ 

-2,949597t 01 

6*03a«35£-02 

B,9CJ0352E 01 

2;965041E^04 

1.576010E 02 

Z;o[i0ia7F .00 

2,296fiflaE 01 

S,Se705B5E-O2 

«5,0b5Sl9E 01 

2‘,339t.a?E 00 

U7B0inaF 02 

«.0376?9B-01 

1 ,493bq0E 02 

i;32752aE 00 

-l*U2q=Jft5E 02 __ 

a,52a753E'*0l . 

-1.72^737E 02 ... 

S;92707bE 00 

jfjS'Toqqe oa 

1,97271BF 00 

3.241450E 01 

t'a9siS7E-oa 

.»2,iJb93ifrt-»0l. 

9^U4abaF-C2 

1,090769£ 02 

6;i65b0bE«03 

5,1&5339E 01 

2,l99ai7£-0l 

-7*B9C653E 01 

i;S3B70BE.-04 

3»‘795e9tjE PI , 

3,Bl6531E-02 

l,oai504E C2 

t;75a770E 00 

0,«5066aF 01 

1:o32089E 00 

q,a74flq4E 01 

sCqOetupaE 00 

l,77bbiB£ 02 

K?2249tF-oa 

, I,6b05q2£, 02 

r«}0fa693E OD 

6,3ai«2SF 00 .. . 

U265051F-12 

-2,B97b60£ 00 

. 2;57l23faF-0i 

• l,«bl2Sta .02 

1.99542BF-01 

l.uaSOSTE 02 

B.HBP0B7F-OB 

-e,S4l9ii5E 01 

l,920lB6t-0b 

-UO'JgTBbF 02 

l,5fl9utt0E-0l 

l^fa7299fiE.C2, 

_i;89512bF.lt . 

.__-5,i469£ib£^Qi„_ 

E'B07b02E*.oa 

l.095a5lE 02 

f^322520E-l4 

0. 

a;2iioa7|E-oi 

• l.lfaT'JtBE 02 

!’,3tlB40E OO 

-B,5at274E 01 

B.H90039E-01 

-l,74aas7F. 02 

l,6704bBE^01 

-7jl90253E 01 

j's5qBb9F-01 . 

a*S2fab|lE 01 . 

lU‘??237F^0ti.. 

-l,06y4<7F 53 

S, 1!5BS7E-o2 

fa,735;«b5 01 

8^l76bl9£^01 

"7.4I7493E 01 

_l'S7SI5gEn0fa 

1,27B221E (U 

1,693B1BE^C1 

a,926494£ 01 

B.52*3a3bE-lO 

7,0b67l2£ QJ 

3.al50l«E*02 

-1,7626R6£ 02 

_..a'.S97t»7£*l0 

B, 197276E 01 ..._ 

5.342273F-02 

i?,37t05lE OO 

S^SS^OPE-Ol 

W35l3l«r 02 

1' afa3064F-02 

3,109Bi3E 01 

l,3lbOSaE-01 

l,26Baq3F 02 

«'739509H-03 

-a,fr00Ba2E 01 

9;77S333R-Q1 

-lt516403E 02 

4.405662E-04] 

9,(UBa23E fll 

_ o;d7S903E-01 

B,9nU7bE 00 

i*705i67E-Qa 

UlOOa99E 02 . 

i;t006SflE-0l 

0,iSl03bE 01 

2,509133E-11 

3,65ZbiBE 01 

5.mi5SF-02._ 

9,qi59qbE_01 

3i554«96E-lft 

0, 

3>0SfaO3E-O2 

9,30013SE 01 

0, 

0* 

:9U79flUF»0a 

; _ 5.6710S3E 0i .. 

r fl43936E-02 

-1,1465338 02 . 

a;B3oesoE^03 

a.qibTbbE 01 

l.90S4i8E-02 

-2,447327E 01 

^1;079B12E*C2 . 

Ul93b7bF.Q2 

9.43772bE-03 

B .723A33E 01 

3,25fca0bE-Q£ 

-a.so'iiict 01 

1^6664035-02 

l,a42b5SE 02 

^;60225bE-Ob 

. •f5tbB7b96E 01 

3,i3629b£-03 

l,356b5faE 02 

2,30l03aF-l0 

•1.391077E 02 

l,Bb0947F"03 

2,6923l4E 01 

0;5qQ732E..U . 

. _’-^5,124‘?5lE 01 

3*906b2fl£-0i 

-_1,341526£ 02 

Z*,SOfl6fllE-Q2 

5,09069*5? 01 

0, 

0. 

5,ll6aq7E-02 

1,611572E 02 

0 .... 

9. 

B;i327b2E"02 

-7,7b97?aF 01 

0. 

o« 

3'b5Bl92F,02 

-'1,£9P060F 02 

9,- 

O4 

i:3700U2E-oH 

-b.qiTlflBE 01 

0, 


5*5‘l09l7F-03 . 

02 

0. 

0, 

fl.95071BE-Oi 

-«,9fabiulE 01 

0. 


_ _ 2;BaflBaaE.02 

-7,BaC37h£ 01 

0, 

0, 

5jiUB«7F-02 

t,3795l9£ 01 

Q, 

9| 

. R;!i?7b3EpQ2 

l,2fl5B53£ 01 

0^ 

0. 

3.b56l92F-02 

i,6492y5E 02 

o; 

0, 

^ l,Oba7q2F-Ott 

*J,6BQ730E 02 


Qi 

3,Z37l72E-0fa 

fl,li?F03lE 01 

q! ~ ^ 

9, 

__ Z,90aO7aE-O7 

. . -i.qabSqSE o2 

0 

._..^4 


71. Modal Participation in First Seven Modes at Each Plane, 
Configuration 10. 
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Table VIII, Attentuation of Individual Forward- Propagating Modes in 
Treated Sections, 


CONFIGURATION 4 CONFIGURATION 7 


Mode 

AdB 

Section 1 

Section 2 

Mode 

AdB 

Section 1 

Section 2 

1 

6.2 

11.7 

1 

5.6 

7.9 

2 

9.0 

22.1 

2 

10.8 

23.1 

3 

10.1 

16-4 

3 

8.8 

27.4 

4 

12.0 

20.6 

4 

11.8 

27.5 


Overall Suppression, 20-6 dB Overall Suppression, 28.3 dB 


CONFIGURATION 8 CONFIGURATION 9 


Mode 

Section 

AdB 

1 Section 2 

Mode 

Section 

AdB 

1 Section 2 

1 

7.9 

5.6 

1 

3.4 

6.6 

2 

23.1 

10.8 

2 

6.7 

18.8 

3 

27.4 

8.8 

3 

5.5 

26.3 

4 

27.5 

11.8 

4 

7.3 

26.8 


Overall Suppression, 17.7 dB Overall Suppression, 17-6 dB 


CONFIGURATION 10 
AdB 

Mode Section 1 Section 2 


1 6.6 11.3 

2 10.5 28.7 

3 8.1 21.9 

4 14.0 23.5 

Overall Suppression, 33.5 dB 
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Table IX. Forward- and Bacltwarcl-Travelling Enei'gy at Duct Planes. 



CONFIGURATION 

_4 


CONFIGURATION 

_9 

Plane 

Forward 

Backward 

Plane 

Forward 

Backward 

Number 

Energy Flux Energy Flux 

Ntimber 

Energy Flux Energy Flux 

1 

0,5978 

-3.618E-3 

1 

6.1216E-3 

-1.7302E-3 

2 

. 0.5978 

-3.618E-3 

2 

6.1216E-3 

-1.7302E-3 

3 

0.5928 

-4.529E-5 

3 

5.9902E-3 

-l,8433E-4 

4 

0.1483 

-5.361E-5 

4 

1.2684E-3 

-3.1554E-4 

:: 5 

0.1522 

-5.508E-6 

5 

1.2931E-3 

-1.5061E-6 

6 

5.136F>3 , 

-3.504E-6 

6 

7.7558E-5 

-4.5399E-6 

7 

5.U2E-3 

0.0 

7 

7.6962E-5 

0.0 

8 

5.112E-3 

0.0 

8 

7.6962E-5 

0.0 


FORWARD aPWL 



FORWARD 

APWL 


Section 1 -6.0 dB Section 1 -6.8 

Section 2 -14.7 dB Section 2 -12.2 



CONFIGURATION 7 


C0NFIGU.RATI0N 8 

Plane 

Forward 

Backward 

Plane 

Forward 

Backward 

Number 

Energy Flux 

Energy Flux 

Number 

Energy Fliix 

Energy Flux 

1 

12.496E-2 

-36.79E-5 

1 

87.3618E-3 

-38.4946E-5 

2 

12.496E-2 

-36.79E-5 

2 

87.3618E-3 

-38.4946E-5 

3 

12.4614E-2 

-10.806E-7 

3 

87.0438E-3 

-54.4828E-8 

4 

17.4547E-3 

-49.9589E-6 

4 

12.5082E-3 

-89.6023E-7 

5 

17.4175E-3 

-23.1133E-9 

5 

12.4874E-3 

-24.7554E-8 

6 

18.4635E-5 

-10.5437E-7 

6 

14.8182E-4 

-11.1900E-6 

7 

18.4154E-5 

0.0 

7 

14.1931E-4 

0.0 

8 

18.4154E-5 

0.0 

8 

14.1931E-4 

0.0 


FORWARD APWL 


FORWARD O.WL 


Section 1 -8 

.5 


Section 

1 -8.4 


Section. 2 -19 

.7 


Section 

2 -9.4 


CONF IGURATION 10 


Plane Foraard Backward 

Number Energy Flux Energy Flux 


1 

2 

3 

4 

5 

6 

7 

8 


2.6626 

2.6626 

2.6057 

5.8379E-1 

0.5900 

1.1755E-3 

1.1560E-3 

1.1560E-3 


-5.3427E-2 

-5.3427E-2 

-1.3690E-3 

-1.9853E-3 

-2.0164E-6 

-3.6073E-5 

0.0 

0.0 


FORVJARD APliTL 

Section 1 -6.6 

Section 2 -27.0 
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At the source of Configuration 7, modes 1, 2, and 4 dominate evenly* 

The energy redistributes into inodes 3 and 4 at Plane 3^ with 3 dominating* 

At the end of the first section, modes 1 and 3 dominate evenly* The energy 
is dominated by mode 3 at Plane 5, but mode 1 dominates by the end of the 
second section. 

Configuration 8 starts with modes 1 and 2 about even at the source. 

The energy redistributes into the first four mode;; at Plane 3, i;ith 1 and 3 
dominating* At Plane 4, only mode 1 retains significant energy* The energy 
redistributes Into modes 1 and 3 entering the second section, and is left in 
modes 1 and 3 at Plane 6, with mode 1 dominating* 

For Configuration 9, the energy starts out in modes 1, 2, and 4, with 
4 dominant. It redistributes into modes 3 and 4 at Plane 3* By Plane 4, 
the energy is in modes 1, 3, and 4* It then redistrihutes into modes 2 and 
3 at Plane 5, with 3 dominant- At the end of the second section, mode 1 
dominates. 

Modes Ij 2j and 4 dominate evenlv at the source of Configuration 10, 
These redistribute into modes 1, 2, and 3 at the entrance to Section ]* At 
the end of the first section, the energy is in modes 1, 2 , and 3, with 1 and 
3 dominating evenly. At Plane 5, the enei'gy goes into modes 3 and 4, witli 
most in mode 3. At the end of the second section, modes 1 and 3 dominate. 

For all configurations except Configuration 9, the effect of reflected 
waves on the overall attenuation is small. In Configuration 9, there is 
appreciable reflected energy at the source plane, accounting for about a 
1,5 dB drop in the effective (measured) suppression. 

The suppression characteristics in the f^lrst section are remarkably 
similar for all configurations, varying from a low of -6.0 dB to a high of 
-8,6 dB in APIJL. Ti’jo points worth noting: Configuration 10, vfnich gives the 
most overall APWL, provides only -6.6 dB in the first section, and the- first 
section of Configuration 8 (which gave “19.7 dB in its reversed position in 
Configuration 7) provides only -8.4 dB v;hen It leads. 

Comparing the ^'least effective suppressor", Configuration 8, to the 
rest, it is apparent that the performance is controlled by the dominant 
participation of the lowest order mode. The redistribution of modes coming 
off the hardwall section is not as advantageous to suppression at it was 
when coming off the "set-up" treatment section in the arrangement of 
Configuration 7. The second section of Configuration 8 is then ineffective 
in redistributing and suppressing the pattern coming from the first section. 

Differences among Configurations 4, 7, 9, and 10 can be explained by 
considering individual mode attenuation rates in the second sections* In 
the second section of Configuration 4, mode 1 has the highest attenuation 
rate of all sections, but at the expense of mode 3, which dominates. In 
Configuration 9, mode 3 has a high .attenuation rate in the second section, 
but mode 1 has a low attenuation rate. Thus, even though the third mode 
dominates at the entrance to Section 2, the suppression is controlled by 



mode 1. For Configuration 7, the overall suppression is aided by a high 
APIIL in the first section. In the second section of Configuration 7, the 
situation is similar to. Configuration 9 except that the mode 1 rate is 
slightly higher. In Configuration 10^ the high suppression in the second 
section appears to be caused by obtaining just the right tradeoff in 
attenuation rates of the first and third modes. At the end of the second 
section 3 in this case* itiodes 1 and 3 are left with approximately equal 
amplitude, which approaches the optimum situation* 

From the above, it appears that the action of the first treated 
segment is to ’*set-up'* a particular pressure profile which the second 
section can effectively "v7xpe-out" at about 4 times the suppression rate 
of the first segment. A possible indication of this trends is apparent 
from the amplitude plots of the pressure profiles at Plane 5. For the two 
high suppression configurations, 7 and 10, there is a pronounced skew in 
level toward the outer wall of the duct. The skew tox^ard the outer x^all is 
less pronounced but nevertheless present in Configurations 4 and 9. In the 
worst performer. Configuration d, the outer wall skex-J is very slight, such 
that most of the energy in concentrated at the center of tue duct . This 
suggests a physical explanation of the operation of a tx;o-element suppressor 
is that the first section distributes the energy flux to the outer x^alls of 
the duct, and the second section tak::±.s advantage of this, situation, as 
discussed in Reference 8. 


4.4.5 Evaluation of Test Results and Correlation. 

Table X summarizes the test configurations for vjhich modal measurements 
x^ere taken and compares the measured and predicted results using nominal imped- 
ance values* The correlation is good for Configurations 4 and 9, but over— 
predicts for Configurations 7, 8, and 10* 

The inability to come x^ithin 10 dB of the predicted optimum attenuation is 
disappointing, and can be traced to tx*7o main problems. First, the modal content 
in the source signal appears to be dependent upon the treatment present and is 
quite sensitive to slight changes in frequency. Second, the predicted atten- 
uation is very sensitive to changes in impedance, particularly near the optimum 
condition. 

The dependence of the source modal content on the treated section could 
possibly be explained by changes in modal pattern at the source plane due to 
modes reflected from the leading edge of the treatment Interacting xotii for- . • 
x^ard-traveling modes. The analytical results, hox^ever, indicate that the 
haclward-traveiing energy ^.t the source plane, at these frequencies, is 

typically 25 dB or more belox-/ the forx>rard-traveiing energy (P\^), indicating 
that their effec.ts should be negligible. It is possible that the. acoustic 
’^loading" presented by the duct upon the source may cause the radiation 
efficiency of different modes to vary from configuration to configuration* 

The sensitivity of the suppression to the changes in modal content x^hich 
Xirere experienced for the different configurations x^as investigated by calcu- 
lating the suppression for a number of different measured sources holding the 
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Table X, Surama^cy of Measured and Predicted Suppressions. 


Cotifigo 



Liner Characteristics 

Impedances 

Z/Pc 

4 

Section 

1 

10% Porosity, 

1.3cra 

deep 

.90 ~1.4i 


Section 

2 

10% Porosity, 

2.3cm 

deep 

.90 -1.7i 

7 

Section 1 

14% Porosity, 

lo8cm 

deep 

.64 -l.lSi 


Section 

2 

12% Porosity, 

2.5ctn 

deep 

.75 -0.55i 

8 

Section 

1 

12% Porosity, 

2 , 5cm 

deep 

.75 -0.55i 


Section 

2 

l47o Porosity, 

1.8cm 

deep 

.64 -1.15i 

9 

Section 

1 

20% Porosity, 

1.8cm 

deep 

.43 -1.3i 


Section 

2 

14% Porosity, 

2.5cm 

deep 

.64 -0,6i 

10 

Section 

1 

12% Porosity, 

1.8cm 

deep 

.75 *rl.2i 


Section 

2 

10% Porosity, 

2. 5cm 

deep 

.90 -0.55i 


Measured 
APWL, dB 




-20o0 


-11. C 


-16.5 


-26.1,-26.2 


Predicted 
APWL, dB 


- 20.6 


-28.3 


-17.7 


-17.6 


-33.5 



impedance components of the treatment constant at the theoretical optimum 
values for the modal content of Configuration 4. Table XI lists the predicted 
suppressions for these cases^ all calculated at 2000 Hz* ]?or these samples j 
any^diere from 2*6 dB to 12*7 dB would be lost due to changes in modal content 
alone if Z/pc were not changed* 

The optimum value of suppression occurs at the top of a very sharp peak 
in the four parameter impedance space of the two liners. This is indicated 
in Table V which shows that variations in any one of the impedance components 
values by 0.1 can cause the suppression to fall off anywhere from 1.2 to 8.4 dB 
from the peak value. Combinations of variations from the optimum impedance 
components could be expected to give further reductions from the peak value. 

Based on the results of the sensitivity studies, it can be concluded that 
the maximum value of suppression practically obtainable under the given two- 
element panel design constraints has been closely approached. The “sharpness" 
of the attenuation peak would seem to indicate that it will be of greater 
' practical value to investigate effects of broadband frequency suppression than 
to attempt to refine the pure- tone optimum design. 



Table XI. Predicted Suppressions for Varying Modal Input, Holding 
Impedance Constant at Optimum Value, Macli 0.3, 2000 Hz, 


MODES OBTAINED FROM 
CONFIGURATION; 

4 (Figui’e 37) 

7 (Figure 40) 

8 (Figure 42) 

9 (Figure 44) 

10 (Figure 46) 


PREDICTED aPWL 
dB 


-35.9 


"27.9 


-23.2 

-28.2 


-23„0 


In all cases 

Z^/Pc = „6 -1.2i 
Z 2 /PC = .7 -.51 

loa 



SECTION 5.0 


ANALYTICAL PARAMETRIC STUDIES 


Results of systematic optimization studies of uniform treatment and of 
two-segment treatment were previously discussed only for a single value of 
n and a single L/H value. In the above two-segment liner configuration the 
segments were of equal length, i.e., - L 2 . It was deemed desirable to 

carry out additional parametric studies. These analytical optimization 
studies included: 

• The effect of r\ (ratio of duct height- to-wavelengtli, H/X) on the 
optimum sound suppression in an acoustically lined rectangular duct 

• The effect of variation of liner segment length ratio on optimum 
sound suppression when the total length of the ti^jo-element treat- 
ment is kept constant. 

The following is a discussion of the results of these parametric studies 


5 . 1 The Effect of n Variation on Optimum Attenuation 

The effect of variation of n on the maximum suppression was explored for 
the following conditions: 

• Rectangular duct geometry modeled in two dimensions. 

• Duct height-to-x^avelength ratios, ri, ranging from 0.5 to 10. 

• Constant total wall liner length with a corresponding L/H of 3.43. 

• In the two-segment liner configuration, a 1:1 ratio of the panel 
lengths was maintained. 

• Exhaust conditions, i.e., flow and sound propagating in the same 
direction. 

• Uniform profile flow of Mach 0,3. 

• Continuity of particle displacement wall boundary conditions. 

• Plane wave source modal input* 

• Three inodes used in calculation* 

The iterative optimization procedure was described in Section 3.2. 
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Table XII summarizes the results of the optimization of uniform treat- 
ment and of two-segment treatment for a range of r\ values from 0.5 to 10. In 
that range of n> optimum sound suppression was predicted for negative . reac- 
tance values. In all cases, the criteria for judging the validity of the 
computer runs was the balance of energy flux at section interfaces. 

Figure 73 graphically depicts the peak suppression and associated opti- 
mum resistance and reactance ratios as a function of n for the case of 
single-element treatment. The negative optimum reactance ratio decreases 
monotonically with increasing ri, i.e. , with increasing frequency. It 
decreases from about -0.1, at n of 0.3, to -3.6 at an q of 10. Tiie. corres- 
ponding resistance ratio increases moderately with n. Optimum suppression 
of the uniform liner decreases markedly with increasing n* 

Optimum suppression predicted for the two-segment treatment and the 
associated liner segment impedance components are shown in Figure 74 as a 
function of n* The optimum resistance of the upstream liner segment is 
relatively small. Its vaJ/ue ranges from about 0.6 pc shown to 0.005 pc for 
high values of q. The optimum resistance of the doi^mstream (second) liner 
increases with Ua i.e., it increases with frequency. The negative reactance 
of both liner segments increase in absolute value with n. Tlie absolute value 
of the reactance of the first segment increases at a faster rate than that 
of the second segment. For X] larger than 1, the negative reactance of . the 
first liner segment is larger, in absolute value, than that of the second 
liner segment. This suggests a smaller cavity depth for the first than that 
of the second optimized liner segment. 

Tlie two-phase treatment optimum suppression also decreases with increase 
in n. Its decrease, however, is not as mono tonic as that of the uniform 
treatment. The too-phase treatment suppi*ession and the uniform treatment 
suppression are compared graphically in Figure. 75. For all considered n 
values the two-phase treatment suppression exceeds that of the single-phase 
liner* The largest beneficial phasing effect of about 10 dB is shown for an 
n of about 1.5, and about 1.0 dB for n values lower than 1.0. 

At then values of 0.5 and 1.0 there is little difference in impedance . 
components bet^^/een the optimized single-element and two-element configura- 
tions. The optimisation becomes extremely difficult in: this x\ region due to 
the high sensitivity of calculated suppression to small changes (on the order 
of O.Ol pc in the. impedance. From practical considerations, it. appears doubt- 
ful that liners could be manufactured to the tolerances required to take 
advantage of the multiphasing effect for these low n values, Gases at these 
low n values were particularly prone to energy flux unbalances, requiring many 
runs, to be discarded. The cause of. these mismatches requires further 
investigation* 

( 12 ) 

Lester and Posey have determined an optimized suppression curve for 
a tx^o-e lament, equal -length liner configuration in a cylindrical geometry 
duct. These curves exliibit a ”knee" phenomenon in the too-element liner 
configuration similar to that of Figure 75, such that the major multiphase 
enhancement occurs in the midrange of n values. 


$ 
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Table XII, Sunliiiary of Optimization Results for Several T| Values, 



w ^ M & SOUND 

Plane Wave Source 


Uniform Liner Two- Segment Liner 


f, Hz 

’I 

Zj^/Pc^ Z^/Pc 

-AdB' 

Z^/Pc 

Z ^ /p c 

-AdB 

650 

0.5 

0.31 -O.li 

70.6 

.32 -0.09i 

0,3 -O.lOi 

71.4 

1303 

1.0 

0.65 -0.5i 

39.2 

.60 -0.47i 

0. 72-0.501 

41.4 

2000 

1.535 

0.85 -l.Oi 

22.6 

.67 -1.37i 

0.60-0.621 

34.3 

2606 

2.0 

0.90 -1.4i 

14.3 

.50 -1.60i 

1.20-0.901 

22.9 

6514 

5,0 

0.85 -2.45i 

4,0 

.005-2.90i 

1.00-1.751 

6.1 

13029 

10 ,0 

1.40 -3.6i 

2.3 

.05 -4.40i 

2,00-3.001 

2,8 


* Z]^/Pc 

:= impedance 
= impedance 

ratio of 
ratio of 

upstream segment 
downs tream s egment 



** "AdB 

= suppression 

in dB's 







Suppress xon 












F i ^ II r e 75 . Max Imiiiu Precl i c t ed S ouncl o i‘ S u ppr s .s ion . i n a Recta n^u 1 a i' Due t w i t h 

■Uni to nil and Two^ Phase Treatment, Mach n, 3, 





5.2 EFFECT QF SEGMENT LENGTH RATIO VAItlATION 

The conditions in this part of the parametric optimization investigation 
were the same as those in the study of the n effect, except as noted below: 

• Duct height-to-wavelength ratio, ri was 1.535. 

« The ratio of the upstream liner segment length (L^) to that of the 
downstream segment (L2) was varied over a range of values from 
1:7 to 2:1, holding the overall length, Lj + L2, constant. The 
duct geometry is depicted schematically in Figure 76. 

0 The source modal distribution was based on measured modal, infor- 
mation from Configuration 4 at Mach 0.3. The modal expansion 
coefficients used here (amplitude and phase) are shown in Figure 77. 

0 The five lowest transverse modes were used. 

The effect of variation of the ratio of the lengths nf the too liner 
segments, with the total liner length kept unchanged, is summarized in . 

Table XIII and shown graphically in Figures 78 and 79. In addition to the 
optimized total suppression. Table XIII and Figure 73 show the associated 
suppressions due to each of the too liner segments. The respective optimized 
resistance and reactance ratios (impedance components) are also included in 
Table XIII and are depicted graphically in Figure 79 . Figures 78 and 79 
show the respective calculated poMt values with the smooth curves drawn 
through them. The uniform (single— phase) liner information, included in Table 
XIII and Figures 78 and 79 were taken from Table II. They were obtained, for 
a slightly different source modal input and seven transverse modes.. 

Table XIII and Figure 78 shot; an increase in total suppression with 
increase in the liner segment length ratio from zero to sli^tly over one. 

The former corresponds to a uniform (one-segment) liner. The length ratio 
of one corresponds to equal lengths of the liner segments. Total optimum 
suppression undergoes little change for the range of segment length ratios 
between 1.6 and 2.0. It appears to have reached its maximum ('v- 36.0 dB) 
at the. segment length ratio of about 1.0, i.e., for the too liner segments 
of equal lengths, as Table XI IT Indicates. Figure 78 shows that it is a 
broad suppression maximum and it extends from the point of length ratio 
(L1/L2) of 1.0 to about 1.25, at which the too liner segments have approxi- 
mately the same ..resistance, , For liner segment length ratio . ahpve about 1.25, 
total suppression decreases with the increase in the segment length ratio. 

Within the considered range, up to the L1/L2 ratio of 2.0, the downstream 
liner se^ent .suppression far exceeds that due to the. upstream segment ^ 

Even at the highest considered ratio of 2, i.e., when the upstream liner 
segment is toice as long as the downstream segment, over toice as much 
predicted suppression is due to the latter ('i^ 24 dB) as that to the former 
('^^ 10 dB) . The suppression due to the upstream liner segment increases very 
moderately with increase in its length at the expense of the length of the 
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f 2000 Hz 
M = 0.3 

L 91c44cm (36 In.) 

H M 26o67ctn (10.5 In.) 

H/X ^ 1«535 


Figure 76. Schematic View of the. Rectangular Duct in Studies of Segment 
Length Ratio Eltect, 



Keiatx' 



MODAL INDEX 


Figure 77, Source Modal Coefficients used in the Studxes of i»iner Segment 
Length Ratio Effect. 




Table XIII. Calculated Effect of Linar L"|ment Length Ratio Variation, 




Z^/Pc 

Z 2 /PC 


-£:.dB2 

-AdB 

0* 

.90 -l.Oi 

.90 -l.Oi 

- 

- 

21.6 

1/7 

.04 -0.9i 

.75 -0.87i 

0.2 

28.0 

28.2 

1/3 

.17 -1.17i 

.75 -0.77i 

1.2 

29.9 

31.1 

1/2 

.30 -1.20i 

.75 -0.68i 

2,8 

29.6 

32.4 

1/1 

.60 -1.20i 

.70 -0.5i 

7.0 

28.9 

35.9 

itx 

.77 -1.17i 

.60 -0.42i 

10.4 

23.9 

34.3 


* Uniform Liner 
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RECmNGULAR DUCT 


f = 2000 Hz 
M = 0.3 

L = Lj^+L^= 36” (91,44 cm) 

HA = 1.535 
L/H « 3.43 


T^tal Suppression. 


•Due to Uniform 
Liner 


Due to Downstream 
Segment 


Due to Upstream 
Segment 


Segment Length Ratio j 


Figure 78. Predicted Effect of Liner Segment Length Ratio on Optimum 
Sound Suppression. 







Reactance Ratio, X/pc Resistance Ratio, R/Pc 






downstream segment. In contrasty the suppression due to the do^^mstream liner 
segment increases rapidly as the single-phase liner is replaced with the ti*;o- 
element treatment, even though the first (upstream) element constitutes only 
one-seventh of the length of the downstream liner segment. That increase 
continues until the upstream segment length is about 25% of the total treat- 
ment length, .i*e. , is about one-third as long as the downstream liner panel. 

The optimised liner segment resistance and reactance ratios are shouTi 
graphically in Figure 79. The optimized reactances of both liner segments 
are negative for the considered segment length ratios- The optimum resistance 
of the upstream segment increases with its length, i.e., with the increase of 
the hx/L 2 ratio. The optimum resistance of the doxmstream segment decreases 
very moderately with the decrease in its length (L 2) 3 i.e., with the increase 
in the Lx/L 2 ratio. The two resistance curves .intersect at an L]^/L 2 ratio of 
about 1^25. As already mentioned, this resistance value appears to be quite 
close to the maximum suppression of the t\^o-element treatment at the considered 
conditions . 

The optimum reactance ratio of the upstream segment is larger, in absolute 
value, than the corresponding reactance ratio of the doximstream panel. This 
indicates a smaller core depth for the upstream liner segment than that for the 
downstream segment in an optimized too-phase liner configuration. 



SECTION 6.0 


CONCLUSIONS 


Aa analysis for calculating the propagation of acoustic waves in rectan- 
gulai- flow ducts has been developed to be used in conjunction with in-duct 
modal measurements. The method has been used to design and test a t\vTo- element 
optimized treatment configuration. 

It has been shown that duct acoustic modes can be measured successfully 
in engine etlvirohment flow and sound conditions , and used to accurately pre- 
dict sound attenuation for various treatment configurations in the presence 
of many modes. Of the five treatment configurations tested, the prediction 
was quite close for two cases, but an undetermined source of error caused 
overpred let ions in the other three. 

The discrete frequency optimum liner designs were found to possess very 
narrow suppression peaks, being quite sensitive to variations in impedance 
components or modal content. This is true for both single-element and two- 
element optimum treatment designs , The sensitivity to impedance components 
and modal content implies a strong dependence of peak suppression on frequency 
variations, as well. At the design Mach number, the measured third-octave 
suppression data for Configurations 4 and 10, however, indicated an enhance- 
ment of broadband suppression for the two-element liner as compared with, 
suppression bandwidth curves for single-element liners obtained by normal 
design procedures, but whether the enhancement was great enough to warrant 
using the two-element design for suppression at a single frequency is 
questionable. 

The highest measured narrowband suppression came within 10 dB of the pre- 
dicted optimum peak suppression of about 36 dB. The combined effects of 
changing source modal characteristics and the inability to precisely achieve 
optimum impedance components accounted for this difference. The two effects 
can be considered to contribute a roughly equal amount to the loss in peak 
suppression. The results of the sensitivity studies indicate that the pre- 
diction technique correlates well with the measured suppression \o.thin the 
uncertainty caused by expected errors in the impedance components* Thus, the 
analysis has been shown to be adequate for the purposes of analytical parametric 
studies . 

The narrowband attenuation obtained for Configuration 10 is likely to be 
as close as the predicted discrete frequency optimum can be approached in 
practical cases. It will be quite difficult to achieve the predicted high 
suppressions within 5 to 10 dB due to the ^’sharpness" of the attenuation peaks 
near optimum, particularly at lower n values. In future studies, it x^ill 
prove more useful to develop this approach for broadband suppression rather 
than refine the discrete frequency development* 
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The results of the analytical study indicate that the two-element liner 
offers enhanced suppression over a single-element liner with the same total 
length. The most practical use of this enhancement occurs over an approximate 
range of frequency parameter q. (H/A) from 1 to 5, For the considered condi- 
tions, the maximum increase in suppression for the two— segment treatment com- 
pared to the single— element treatment was predicted to occur at an p value of 
about 1.5, Analysis of the computer runs indicates that the primary two- 
phase mechanism causing the enhanced suppression at this p value is modal 
redistribution from one plane to the next. The effect of energy reflection at 
the interfaces is negligible. 

Comparison of the experimental results for Configurations 7 and 8, with 
the same panels in reversed positions, indicate that the two liners do work in 
conjunction with one another. The design of the upstream treatment section is 
dependent on the design of the second; that is, one could not optimize the design 
of the first section for a given modal input, then use the output of the first 
section to optimize the design of the second section and expect to get the 
maximum suppression of both sections taken together. The optimization studies 
have Indicated that several sequential iterations on both sections are 
necessary before the overall suppression converges to the peak value. The 
studies have indicated that, holding all other parameters except the tt^o 
liner impedances constant, there is a unique set of the four impedance com- 
ponents which maximizes overall suppression. 

Results of the analytical parametric studies have shown that maximum sound 
attentuation due to a two-segment liner of a given length depends upon the rela- 
tive length of the liner segments. Optimum suppression increases with increase 
from zero to slightly over one in the ratio of the length of the upstream 
liner segment to that of the doximstream segment. Maximum optimized suppres- 
sion is attained at or slightly above the segment length ratio of 1.0, i.e,, 
when the tii?o liner segments are approximately of the same length. For the 
considered conditions, up to the L]^/L 2 ratio of 2, the downstream liner seg- 
ment suppression by far exceeded the suppression due to the upstream segment. 

The modal content of the duct source varies rapidly xfith frequency for 
pure tone excitation. Further investigation is needed to determine whether 
the change of the source modal content x-irith frequency and in the presence of 
different treatment configurations xjas just a characteristic of the rectangu- 
lar duct apparatus or xv’hether it Is of significance in engine ducts as xvell. 

The analytical prediction procedure using measured modal content has 
been shown to be an adequate design device for discrete frequencies. In 
future studies, it is recommended that the method be adapted for predictions 
over a broadband range of frequencies. This would be of particular importance 
in investigating such effects as the enhanced high frequency too- element liner 
suppression measured in this study. It xjould be useful for the development 
of broadband treatment design criteria. 
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SECTION 8.0 


LIST OF SYMBOLS 


+ 

A" 

m 



ib 

E (mjn) 
z 


= modal expansion coefficients 

= modal solution vector at plane j 
= speed of sound in stationary medium 

= a"!” + a . 

= modal energy flux 


f = frequency, Hz 

F = general function 

H = duct height 

i = /=! 

- the axial component of the acoustic intensity vector 
k = the wave number in free space 

I - liner cavity depth 

M = axial floxir Mach number 

n = index indicating type of boundary condition 

NM = number of transverse modes in the expansion 

p(x,Zjt) “ the acoustic pressure at point (x^z) and time t 
{Q} = the stacked modal source vector 

R/pc = resistance ratio 

^ reflection matrix at plane j + 1 of the forward 
traveling wave 

SCoJax) " cross- spectral density 

[S] - the system matrix 


t = time 

tp = face plate thickness 

_ transmission matrix from plane j to plane j -f 1 for 
the fo'n'yard propagating wave 

U = mean flow velocity in the axial direction 

[U] ” uniform section transmission matrix 


122 



LIST OF SYMBOLS (Concluded) 


Xj z 

X/pc 

Z/pc 

3 

TTm 

6 

^IJ 

n = H/A 

+ 

K" 

m 

X 

P.PQ 


(Ai 


- axial component of acoustic velocity 

- Carthesian Coordinates in too-dimensions 

- reactance ratio 
” impedance ratio 
= admittance ratio 

= eigenvalue of the M~th transverse mode 

- empirical liner hole end correction 

- the Kronecker Delta 

- duct height-tO“Wavelength ratio 

= axial propagation constant for the M-th mode 

= wavelength 
= density of air 

- faceplate open area ratio 

- real part of propagation constant, m*^h mode 

- imaginary part of propagation constant, mode 

" M-th mode eigenfunction 

= vector of eigenfunction at plane j 
_ 2Trf = circular frequency 
= the Laplacian Operator 


m. . 

P 

r 

T 

i 

* 

< > 



axial station indices 
transverse mode number 
referring ho probe 
referring to reference microphone 
designating a transpose matrix 

superscript sign denoting the direction of wave 
propagation 

designates complex conjugate 
brackets designating time average 
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